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Systematic  studies  of  the  nuclear  magnetism  of  solid  ortho-para  hydrogen 
mixtures  at  low  temperatures  are  presented. 

The  formulation  of  the  nuclear  spin-lattice  relaxation  time  Ti  for  the  case 
of  the  local  ordering  in  the  orientational  "glass"  phase  of  ortho-para  hydrogen 
mixtures  is  given.  The  temperature  dependence  of  Ti  is  discussed.  A  strong 
dependence  on  the  position  of  the  NMR  isochromat  in  the  line  shape  is  found 
and  this  is  in  good  agreement  with  the  experimental  results  of  a  group  of 
physicists  at  Duke  University,  provided  that  the  cross-relaxation  is  taken  into 
account.  The  relaxation  is  found  to  deviate  considerably  from  an  exponential 
recovery. 

The  orientational  degrees  of  freedom  of  ortho  hydrogen  molecules  in  terms 
of  density  matrices  and  the  irreducible  tensorial  operators  associated  with  unit 
angular  momentum  are  described.  The  range  of  allowed  values  for  the  orienta- 
tional order  parameters  is  determined  from  the  positivity  conditions  imposed 
on  the  density  matrix. 

A  theory  of  solid  echoes  and  nuclear  spin  stimulated  echoes  following  a  two- 
pulse  RF  sequence  and  a  three-pulse  sequence  in  the  quadrupolar  glass  phase 
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of  solid  hydrogen  is  developd.  The  stimulated  echoes  can  be  used  to  compare  a 
"fingerprint"  of  the  local  molecular  orientations  at  a  given  time  with  those  at 
some  later  time  (less  than  Ti  )  and  thereby  used  to  detect  ultra-slow  molecular 
re-orientations. 
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CHAPTER  1 

INTRODUCTION 

The  orientational  ordering  of  the  rotational  degrees  of  freedom  of  hydrogen 
molecules  at  low  temperatures  has  been  carried  out  by  a  number  of  groups  in 
recent  years  (1  -  19).  The  principal  reason  for  this  interest  is  that  the  ortho 
hydrogen  (or  para  deuterium)  molecules  with  unit  angular  momentum,  J  =  1, 
represent  an  almost  ideal  example  of  interacting  "spin-l"  quantum  rotators, 
and  therefore  a  valuable  testing  ground  for  theoretical  models  of  co-operative 
behavior. 

A  suggestion  that  a  random  distribution  of  ortho  molecules  in  a  para- 
hydrogen  matrix  may  behave  as  a  quadrupolar  glass  at  low  temperatures  for 
ortho  concentration  X  <  55%  was  proposed.  It  was  based  on  the  observa- 
tion that  ortho-para  hydrogen  mixtures  provide  a  striking  physical  realization 
of  the  combined  effects  of  frustration  and  disorder  on  collective  phenomena. 
These  effects  play  a  determining  role  in  spin  glasses,  and  the  behavior  of  solid 
hydrogen  mixtures  is  analogous  to  that  of  a  spin-glass  such  as  EuxSri_^S  in 
a  random  field(20,21). 

The  electrostatic  quadrupole-quadrupole  interaction  is  the  dominant  inter- 
action which  determines  the  relative  orientations  of  the  ortho  molecules,  and 
there  is  a  fundamental  topological  incompatibility  between  the  configuration 
for  the  lowest  energy  for  a  pair  of  molecules  (a  Tee  Configuration  for  EQQ) 
and  the  crystal  lattice  structure:  one  cannot  arrange  all  molecules  so  that  they 
are  mutually  perpendicular  on  any  3D  lattice.  '■ 
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At  high  temperatures  the  molecules  are  free  to  rotate,  but  on  cooling, 
the  ortho  molecules  tend  to  orient  preferentially  along  local  axes  to  mini- 
mize their  anisotropic  EQQ  interactions.  As  the  temperature  is  reduced,  this 
leads  to  a  continuous  but  relatively  rapid  growth  of  local  order  parameters, 
[a  =  (3J^  -  2))  which  measure  the  degree  of  alignment  along  the  local  sym- 
metry axes  oz.  There  is  a  broad  distribution  p(a)  of  local  order  parameters  at 
low  temperatures  (12,22,23),  but  no  clear  phase  transition  has  been  detected 
(18,24). 

The  degree  of  cooperativity  in  the  slowing  down  of  the  orientational  fluctua- 
tions of  the  molecules  as  the  samples  are  cooled  is  of  special  interest  in  these  sys- 
tems and  this  has  motivated  several  independent  experiemental  studies(25-28) 
of  NMR  relaxation  times  Ti,  which  are  determined  by  the  fluctuations  of  the 
molecular  orientations.  One  of  the  most  striking  results  reported  by  S.  Wash- 
burn et  al.  (28)  is  the  observation  of  a  very  strong  dependence  of  Ti  on  the 
spectral  position  within  the  NMR  absorption  spectrum.  While  a  qualitative  in- 
terpretation of  this  behavior  in  terms  of  "sloppy"  librons  has  been  off"ered(l,28), 
a  more  detailed  treatment  has  been  lacking.  One  of  the  aims  of  my  research  is 
to  extend  earlier  work  (29)  by  using  a  straightforward  theory  and  to  compare 
the  results  with  the  experiemental  data. 

In  the  low  concentration  regime,  although  it  is  agreed  that  (i)  the  low 
temperature  NMR  lineshapes  indicate  a  random  distribution  of  molecular  ori- 
entations (for  both  the  local  axes  and  and  the  alignment  a  =<  3j|  —  2  >,  and 
(ii)  that  there  is  apparently  no  abrupt  transition  in  the  thermodynamic  sense; 
there  has  been  disagreement  (1,3)  over  the  interpretation  of  the  behavior  of 
the  molecular  orientational  fluctuations  on  cooling  from  the  high  temperature 
(free  rotator)  phase  to  very  low  temperatures  (T  <  O.IK).   While  some  early 
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work  reported  a  very  rapid,  but  smooth  variation  (25,27,30)  with  tempera- 
ture, corresponding  to  a  collective  freezing  of  the  orientational  fluctuations, 
subsequent  studies  (i4, 26,28)  indicated  a  slow,  smooth  dependence  with  no 
evidence  of  any  strong  collective  behavior.  In  order  to  resolve  this  problem  it 
is  important  to  understand  two  unusual  properties  of  the  nuclear  spin  relax- 
ation rates  in  the  glass  regime.  These  two  properties  which  will  be  discussed 
in  Chapter  3  are  (i)  the  spectral  inhomogeneity  (14,29,30)  of  the  relaxation 
rate  T^  across  the  NMR  absorption  spectrum,  and  (ii)  the  nonexponential 
decay  (14,31)  of  the  magnetization  of  a  given  isochromat;  both  result  from  the 
broad  distribution  of  local  axes  and  alignments  for  the  molecular  orientations. 
Both  properties  alone  lead  to  variations  by  more  than  an  order  of  magtitude 
and  need  to  be  understood  theoretically  before  attempting  to  deduce  charac- 
teristic molecular  fluctuation  rates  from  the  relaxation  times.  It  will  be  shown 
that  the  strong  departure  from  exponential  decay  for  the  magnetization  can 
be  understood  provided  that  the  broad  distribution  of  local  order  parameters 
is  correctly  accounted  for. 

In  the  so-called  "quadrupolar  glass"  the  quantum  rotors  cannot  in  gen- 
eral be  described  by  pure  states  and  a  density  matrix  formalism  is  needed  to 
describe  the  orientational  degrees  of  freedom.  It  is  needed  to  determine  the 
precise  limitations  on  the  local  order  parameters  (molecular  alignement,  etc.) 
from  the  quantum  mechanical  conditions  imposed  on  the  density  matrix  and 
to  discuss  the  implications  for  the  analysis  of  NMR  experiments.  Some  of  the 
considerations  for  the  spin-1  density  matrix  description  have  been  given  else- 
where (29, 32  —  37)  but  solid  H2  is  a  special  case  because  the  orbital  angular 
momentum  is  quenched.  The  case  will  be  discussed  in  chapter  4. 
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In  molecular  solids,  one  has  to  consider  two  kinds  of  degrees  of  freedom: 
translational  degrees  of  freedom  for  the  centre-of-mass  motion  and  orientational 
degrees  of  freedom  for  the  rotational  motion  of  the  molecules. 

One  of  the  most  fascinating  problems  encountered  in  the  molecular  solids 
is  the  existence  of  glass-like  phases  in  which  the  molecular  orientations  become 
frozen  without  any  significant  periodic  correlation  from  one  site  to  another 
throughout  the  crystal.  The  most  striking  example  of  these  "orientational 
glasses"  is  probably  that  observed  in  the  solid  hydrogen  when  the  quadrupole- 
bearing  molecules  are  replaced  by  a  sufficient  number  of  inert  diluants. 

If  the  quadrupoles  in  this  frustrated  system  (20)  are  replaced  by  "inert" 
molecules,  it  leads  to  large  reorientations  of  the  quadrupole-bearing  molecules 
in  the  neighbourhood  of  the  inactive  diluants  and  the  dissappearance  of  long- 
range  order  when  the  quadrupole  concentrations  is  reduced  below  55%.  The 
HCP  lattice  is  apparently  stable  down  to  very  low  temperatures  (7)  and  the 
NMR  experiments  indicate  that  the  molecular  orientations  vary  in  a  random 
fashion  from  one  site  to  another,  both  the  directions  of  the  local  equilibrium 
axes  and  the  degree  of  orientation  with  respect  to  these  axes  vary  randomly 
throughout  a  given  crystal. 

The  important  questions  are  (i)  whether  or  not  the  freezing  of  the  molecular 
motion  persists  for  time  scales  much  longer  than  those  previously  established 
for  the  glass  phase  and  (ii)  how  the  freezing  occurs  on  cooling. 

In  order  to  answer  these  questions  a  new  type  of  experiment  was  clearly 
needed. 

The  echo  techniques  are  of  great  practical  importance  in  NMR  measure- 
ments on  the  orientationally  ordered  hydrogen,  because  they  allow  one  to  ex- 
tract information  which  is  not  easily  or  unambiguously  determinable  from  ei- 
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ther  steady-state  line  shape  or  FID  analyses.  Conventional  continuous  wave 
(CW)  and  free  induction  decay  (FID)  NMR  techniques  have  only  been  able  to 
show  that  the  orientational  degrees  of  freedom  appear  to  be  fixed  for  times  up 
to  lO""*  —  10~^  S.  A  considerable  improvement  may  be  achieved  by  the  analy- 
sis of  solid  echoes  for  which  the  observation  times  during  which  one  can  follow 
molecular  reorientations  are  extended  to  an  effective  relaxation  time  {T2)effi 
which  may  (by  a  suitable  choice  of  pulses)  be  much  longer  than  the  transverse 
relaxation  time  J2  that  limits  the  conventional  techniques. 

It  will  be  shown  that  spin  echoes  and  stimulated  echoes  following  a  two- 
pulse  sequence  and  a  three-pulse  sequence,  respectively,  provide  a  more  pow- 
erful means  of  investigating  the  orientational  states  and  particularly  the  dy- 
namics of  the  molecules  bearing  the  resonant  nuclei,  than  the  conventional 
continuous-wave  technique. 

Spin  echoes  were  observed  in  solid  H2  a  long  time  ago  (38)  and  have  been 
used  to  study  the  problem  of  orientational  ordering  (11,  14,  27,  28,  39,  40).  In 
order  to  gain  a  deeper  insight  into  this  problem,  a  series  of  questions  will  be 
disscused  in  chapter  5.  These  are  the  formation  of  spin  echoes  (including  solid 
echoes  and  stimulated  echoes),  explanation  and  comparison  with  experiments, 
and  the  motional  damping  of  echoes. 


■CHAPTER  2  ; 

THEORY  OF  NMR  RELAXATION 

Bloembergen-Purcell-Pound  Theory 

The  theory  of  spin  relaxation  in  liquids  (or  gases)  is  based  upon  time- 
dependent  perturbation  theory. 

In  liquids  where  the  spin-spin  coupling  is  weak  and  comparable  to  the 
coupling  of  the  spins  with  the  lattice,  it  is  legitimate  to  consider  individual 
spins,  or  at  most  groups  of  spins  inside  a  molecule,  as  separate  systems  coupled 
independently  to  a  thermal  bath,  the  lattice. 

It  is  well  known  that  the  expression  for  the  interaction  between  two  mag- 
netic dipoles  of  nuclear  spin  I^  and  Ij  can  be  expanded  into  6  terms: 

Hdd  =  ^-^^[h  ■  Ij  -  Hir  ■  n){I,  ■  n)] 
^'^'^^^{A  +  B  +  C  +  D  +  E  +  F)  (2.1) 


where 


C  =  -\{ltljz  +  kJJ)sin6cosde-''^ 
D  =  -^(/r^y^  +  IiJ-)sinecosee''t> 

4  *    J 
4  «    J 
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If  we  carry  out  first-order  time-dependent  perturbation  calculations  to  ob- 
tain the  transition  probability  W  between  tiie  megnetic  energy  levels,  we  find 
for  a  nuclear  spin  /  =  j- 

W  =  -^n^hh[I  +  l)[JiK)  +  ^  J2(2t^o)]  (2.2) 

and  the  longitudinal  relaxation  time 


Ti       2  '         '  yi   XV  w   .   2 

This  is  the  BPP  expression(41,42);  where 

'+00 


2W 
l^h'Iil  +  l)\Ji{uo)  +  ^-J2[2uq)\  (2.3) 


1        3.2 1 


/+00 
-oo 


G{r)  =F{t)F*{t  +  T) 

jPq  =     (1  _  scos'^e) 

Fi  =  -^sinOcosOe^''' 
F2  =  ^sin^Be^'^ 

J.6 


Theory  for  Liquids  in  Terms  of  Mori's  Formalism 
Daniel  Kivelson  and  Kenneth  Ogan  reformulated  the  study  of  spin  relax- 
ation in  liquids  in  terms  of  Mori's  statistical  mechanical  theory  of  transport 
phenomena(43,44,45,46).  They  started  with  some  well-known  phenomenologi- 
cal  magnetic  relaxation  relations  and  formulated  them  in  a  manner  most  suit- 
able for  comparison  with  Mori's  theory.  They  obtained  simple  Bloch  equations 
by  the  Mori  method  and  extend  the  treatment  to  a  time  domain  not  adequately 
described  by  the  Bloch  equations.  For  this  dissertation  I  will  just  demonstrate 
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the  theory  for  the  simple  case  in  which  the  Bloch  equations  can  describe  the 
relaxation  phenomena. 

We  know  for  a  Brownian  particle  that  the  Langevin  equation  is  a  valid 
equation  of  motion  for  times  much  longer  than  the  characteristic  molecular 
times: 

m—  =  I-^v  +  F{t)  (2.4) 

Where  v  is  the  velocity  of  the  Brownian  particle,  7  is  the  force  on  it  due  to 
an  externally  applied  field,  fv  represents  the  slow,  frictional  force  where  f  is 
the  friction  constant,  and  F{t)  is  a  force  which  is  a  rapidly  varying  random 
function  which  averages  to  zero. 

Mori  developed  a  "generalized  Langevin  equation"  which  provides  a  divi- 
sion of  time  scales  into  a  slow  time  scale  associated  with  the  motion  of  the 
Brownian  particle  and  a  fast  time  scale  associated  with  collective  motion. 

Mori  chose  a  set  of  dynamical  variables,  for  example,  A{t),  which  describe 
the  relevant  slow  variations  in  the  system.  A{t)  describes  a  displacement  from 
equilibrium,  i.e.  {A{t))  —  0. 

Let  a  represent  component  of  the  time  derivative  of  A  that  are  orthogonal 
to  A,  and  are  rapidly  varying 

'a  =  {\-P)A  (2.5) 

where  P  is  a  projection  operator. 

The  time  dependence  of  A{t)  can  be  expressed  in  terms  of  the  superoper- 
ator,  ^^,  of  the  Hamiltonian  as 

A  =  {){=' A  '  (2.6) 

A{t)  =  e'^''A  (2.7) 
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The  time  evolution  of  a{tp)  is  determined  by  a  propagator  composed  only 
of  those  components  of  the  Hamiltonian  ){^,  which  lie  outside  the  subspace 
determined  by  the  slow  variables,  i.e. 

k{tp)  =  exp[t{l  -  P)i)^'']k  (2.8) 

We  now  define  a  memory  function  matrix  k{t),  an  effective  memory  func- 
tion  matrix  K{t)  and  a  relaxation  matrix  K  as  follows: 
Memory  function  matrix: 

k{t)  ={k{tp)k)-{AA)-'^  '  (2.9) 

Effective  memory  function  matrix: 

K{t)  =  k{t)e-'^'  (2.10) 

where 

in={AA^)-{AA^)-'^  (2.11) 

Relaxation  matrix  K: 

poo 

K  =  /      K{t)dt  (2.12) 

The  relaxation  matrix  may  be  complex: 


K  =  Kre-  iKiM  (2.13) 

The  real  part  is  associated  with  relaxation  times  and  the  imaginary  part  with 
frequency  shifts. 

Substituting  (2.13)  into  generalized  Langevin  equation,  we  obtain  the  fun- 
damental relation  -  Mori's  equation: 

±A{t)  =  -\-i{n  +  Kim)  +  Kre]  •  W)  (2.i4) 
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and  its  Fourier  transform: 

iu}A{oj)-A{0)=in-A{oj)-K{uj-n)-A{uj)  (2.15) 

We  can  now  derive  Bloch's  equation  and  relaxation  time  formulae  in  a 
simple  case  (one  variable) . 

We  express  the  Ilamiltonian  as 

X  =  XS  +  ^L  +  ^SL  (2.16) 

where  ^5  depends  only  on  the  spin  variable,  "Hi  represents  the  molecular  mo- 
tions and  interactions  which  are  independent  of  the  spins,  and  "Hgi  involves 
those  interactions  which  involve  both  spins  and  nuclear  spatial  coordinates. 
The  slow  variables  A  can  be  selected  as: 

A  =      AS^  (2.17) 

We  assume  that  k[i),  the  memory  function  matrix,  decays  rapidly  so  that 
/C(w  -  n)  can  be  replaced  by  Kre  "  i^IM-  Equation  (2.15),  transformed  back 
to  the  time  domain,  then  becomes: 

d 


^^S,{t)  =  -T-\S,{t)-{S,)]  (2.18) 

-S^  =  [i{±uo  +  a±)  -  T2^]S±{t)  (2.19) 


dt 
dt 


where 


Ti 


4     f°° 

-1  =  ReK,,  =  -  J      {[)isL{tp),Sz{tp)\[S^,MsL])dt  (2.20) 

2      f°° 
=  M^)J^     {{XsL{tp),S±{tp)][S^,MsL])e^''''*  dt  (2.21) 

2      f°° 

^±  =  {KlM)±±  =  rm{j^)J^     {[MsL{tp),S±{tp)]\S^,)^SL])^^""''dt      (2.22) 

Equations  (2.18)  and  (2.19)  are  the  simple  Bloch  equations  except  that  the  a± 
term  which  represents  the  so-called  nonsecular  or  dynamic  frequency  shift  is 
given  explicitly. 
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,    .       Kubo  and  Tomita  Theory 

The  theory  originally  introduced  by  Kubo  and  Tomita  emphasizes  the  simi- 
larity of  magnetic  relaxation  to  other  non-equlibrium  phenomena  (46,47,48,49). 
Description 

We  know  that  the  Bloch  equation  can  describe  the  relaxation  of  the  mag- 
netization m.  It  is  a  linear  equation. 

^  =  ^rn{t)  X  H{t)  -  ^-(^)^^+  ^v(0/  _  m.{t)-m0.j^  (2.23) 

In  matrix  notation  the  Bloch  equation  (2.23)  can  be  written  as 

Ak{t)  =  -L  •  Am(«)  (2.24) 

where 

Arrikit)  =  mk{t)  -  m^     (A;  =  i,j,k)  ]  (2.25) 

m^  is  the  thermal  equilibrium  value 

/-^  0  0         \ 

-L=         0        -^  -iuQ  0  (2.26) 

V    0  0  -^  +  iuoJ 

If  we  transform  from  the  laboratory  frame  to  a  frame  rotating  with  the 
Larmor  frequency  around  the  z  axis,  we  will  have 

Ain{t)  =  -L'Am^{t)  (2.27) 

The  formal  solution  of  Eq.  (2.24)  is  the  following: 

Ama{t)  =  (e-^^)^^Am^(f)  (2.28) 

The  Fourier  transform  of  the  formal  solution  is  therefore 

Ama{iv)+  =  (-^-^)«^Am^(0)  (2.29) 


12  i 


+  means  a  positive  Fourier  transform  and  in  the  rotating  frame: 

Amj{uj)+  =  (— _— )Am^(0)  (2.30) 


'  L-    —  iu 
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Response  Function 


The  linear  response  formalism  begins  by  calculating  the  linear  response 
of  a  dynamic  variable  to  a  disturbance  created  by  a  time-dependent  external 
force. 

Consider  the  Hamiltonian: 

Htotalit)  =  H  -  ^MiHiit)  (2.31) 

where 

(i)H  describes  all  the  interactions  responsible  for  the  motions  of  the  spins, 
including  the  effect  of  the  large,  static  Zeeman  field. 

{n)Hi{t)  is  a  small,  time-varying  field  which  is  responsible  for  the  nonequlib- 
rium  behavior  of  the  system. 

The  response  function  fizi{t)  is  defined  as 

Amfc(^)  =   /      fkiit  -  T)Hi[T)dT  (2.32) 

J  —  oo 

The  response  function  fki[t)  gives  the  effect  of  the  disturbance  at  time  t. 
The  Fourier  transform  of  eq.  (2.32)  is 

^rn^ioj)  =  Xkl{^)Hi[uj)  (2.33) 

where 

roo 

Xkl{<^)  =  /      e''''fki{t)dt  (2.34) 

Jo 

Equation  (2.34)  defines  the  susceptibility  Xkli'^)-  I*  has  a  real  and  imaginary 
part  that  are  connceted  by  the  Kramers-Kronig  relations.  The  imaginary  part 
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of  the  susceptibility  is  called  the  absorptive  part  which  is  related  to  the  power 
the  sample  absorbs.  The  real  part  of  the  susceptibility  is  called  the  dispersive 
part  and  is  related  to  the  measured  line  shape. 
Relaxation  Function 

The  relationship  between  the  relaxation  function  Fi^i{t)  and  the  response 
function  fki{t)  is  given  by 

fkl{r)dT  (2.35) 

The  relaxation  function  describes  the  time  change  of  the  response  after  the 
external  disturbance  is  cut  down  to  zero. 
Assume  a  step  disturbance: 

Hi{t)  =  Hie''e{t)  (2.36) 

where 

and  e  is  a  small  positive  constant  that  will  be  taken  to  zero  at  the  end  of 
the  calculation.  In  principle  this  disturbance  corresponds  to  having  a  field  in 
addition  to  the  Zeeman  field. 

The  response  to  the  step  disturbance  in  the  limit  e  — >  0  is: 

^rnkit)  =  Fki{t)Hi     t>0  (2.38) 

and 

Amfc(O)  =  Fki{0)Hi  ,  (2.39) 

Here  Fki{t)  is  the  relaxation  function  which  describes  how  the  response  to  a 
step  function  disturbance  decays  in  time. 
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If  we  regard  equations  (2.38)  and  (2.39)  as  matrix  equations  and  formally 
eliminate  the  external  force  in  these  two  equations,  we  obtain 


Ama{t)  =  Fa^{t)F^^^O)Am^{0)    t  >  0  (2.40) 

The  central  assumption  of  the  linear  response  theory  is  that  Eqs.    (2.40) 
and  (2.28)  can  be  combined  to  yield  a  molecular  expression  for  L  as 

[e-^'U  =  Fa^{t)F;^H0)  (2.41) 

We  define  a^/(w)  as  the  following 

From  equations  (2.38)  and  (2.39): 

Ama(u;)  =  o^p[uj)Hp  (2.43) 

Ama(O)  =  Xap{^)Hp  (2.44) 

Ama(u;)  =aa^(w)x~j(0)Am^(0)  (2.45) 

Comparing  equation  (2.28)  with  equation  (2.45)  yields 


in  the  rotating  frame 


'kk 
i.e. 


^'kk       _  Re(7i,_k{u  +  Kuq) 


L',\  +  ^'  Xk-k{0) 

Using  a  symmetric  form  of  ajt_jt,  equation  (2.48)  becomes 

^'kk     ^  R^<^k~kH 

L'^l^  +  oj^         Xfc-fc(O) 


(2.46) 


77 -  =  .,  k  =  0,±l  2.47 


(2.48) 


(2.49) 
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Time  correlation  Function  Formulas  for  Transport  Coefficients        ' 

The  transport  coefficient  is  expressed  as  a  time  integral  of  a  correlation 
function  of  magnetization.  These  formulas  will  be  derived  in  the  limit  of  weak 
coupling  between  the  spin  and  lattice  degrees  of  freedom. 

We  assume  that  the  Hamiltonian  of  the  system  H{X)  can  be  split  into  two 
parts:  a  part  //q  that  contains  the  Zeeman  Hamiltonian  and  a  part  H'  that 
couples  the  spins  to  the  lattice  and  is  responsible  for  the  relaxation: 

H{X)  =  Ho  +  \H'  (2.50) 

We  also  assume  the  transport  coefficient  L'[X)  may  be  developed  in  a  power 
series  in  A  with  a  leading  term  in  A'^: 

^'('^)  =  A2X]A"L'(fi)  (2.51) 

Since  we  are  assuming  weak  coupling,  we  identify  the  measured  transport 
coefficient  with  the  first  term  of  the  sum  in  Eq.  (2.45),  i.e.  A^L'(O)  =  -L'(A). 
The  results  of  transport  coefficients  are 

^°^  ^  7\  ^  2(^A^ /o     '^^^\H\M,]\H'[t),M,])^  +  C.C.  (2.52) 

1 


-^11  =  ■C'-i-i  =  ;pr 
-'2 


i\2 


ik) 


where 


roo 

+  {[H'{t),M^][H'  ,M+])q  +  C.C.}  (2.53) 


Xr2  =  ^([AM-,AM+]+>o  (2.54) 

H'{t)  =  eh^^^H'e-hMt)  (2.55) 
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The  subscript  zero  on  the  braket  indicates  that  the  trace  is  taken  over  the 
equilibrium  density  matrix 


exp{-f3HQ) 
r{exp{-l3Hc 

which  does  not  involve  the  spin-lattice  coupling. 


""  =  Tr(eM-m))  <'■''' 


Nuclear  Spin-Lattice  Relaxation  in  Non-Metallic  Solids 
The  problem  here  is  essentially  the  same  as  that  for  liquids  and  gases, 
namely  to  calculate  the  probability  of  a  flip  of  a  nuclear  spin  caused  by  its 
coupling  with  the  thermal  motion  of  a  "lattice."  There  are,  however,  some 
significant  differences.  The  internal  motion  in  solids  will  often  have  much 
smaller  amplitudes  and/or  much  longer  correlation  times  than  in  liquids.  In 
rigid  solids  because  of  the  tight  coupling  between  nuclear  spins  exemplified  by 
frequent  flip-fiops  between  neighbours,  the  correct  approach  to  nuclear  mag- 
netism is  a  collective  one,  where  single  large  spin-system  with  many  degrees 
of  freedom  are  to  be  considered,  rather  than  a  collection  of  individual  spins. 
The  assumption  is  usually  made  that  the  strong  coupling  of  the  nuclei  simply 
establishes  a  common  temperature  called  a  "spin  temperature,"  and  that  the 
lattice  coupling  causes  this  temperature  to  change(50,51,52). 
A  quite  general  equation  can  be  derived. 

i  =  -fli»-M  (2.57) 

where 

0  =  j^Y~    -^s  "  ^P''"  temperature 

00  —  j(Y    "^  ~  Temperature  of  the  lattice 
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Under  certain  conditions  (practically  all  experimental  situations),  high  spin 
and  lattice  temperature,  Abragam  and  Slichter  (51,52)  give  the  following  result 

1    _li:n,mWmn{En-Em)^ 

Ti       2  {M^)  ^2.58) 

Where  |  m)  and  |  n)  are  the  eigenstates  of  ^q- 

Wmn  =  Wnm  is  the  transition  probability  from  the  state  |  m)  to  the  state 
I  n). 

There  is  another  way  to  calculate  Ti.  It  is  a  density  matrix  method,  which 
is  quite  general  and  especially  suited  to  discussing  cases  in  which  motional 
narrowing  takes  place.  Let  the  density  matrix  p  describe  the  behavior  of  the 
combined  quantum  mechanical  system,  spins  +  lattice.  In  the  interaction 
representation 

p*  =  e~ft^«*peh^"^  (2.59) 

7^  =  -[^r(0,/]  (2.60) 

Where  ^i  is  a  perturbation, 

Equation(2.60),  integrated  by  successive  approximation,  gives 

~  =  -{irnhp^m  -^j^dt'mt),{)i*{t'),P*m\ 

+higher  order  terms  (2-61) 


or 


dp*  i 

~dt    ~~h 


\m)^p*m  -  A  /' dT\r,[t), [ui{t - r),/(o)]] 

h    Jo 
-\-higher  order  terms  (2.62) 
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Since  all  the  observations  are  performed  on  the  spin  system,  all  the  relevant 
information  is  contained  in  the  reduced  density  matrix  a* 

a*  =  trf{p*}  (2.63) 

with  matrix  elements  (a  |  a*  |  a')  =  ^f{fa  |  p*  |  fa'). 

By  making  some  assumptions  Abragam  (51)  gives  a  general  master  equa- 
tion 

da*  f°° 

where  the  bar  represents  an  average  of  many  particles. 
If  there  is  a  spin  temperature,  then 

dd       1  r+oo 


I.e. 


ft  = -\%2^  JZ  ^^^*^^'>^^^'^^^*^^^  (2.66) 

1       11     r+°° 

?\^  2(^)7-00    W(0,^o][^r(«-O,>/ol)rfr  (2.67) 

For  the  relaxation  of  like  spins  by  dipolar  coupling  the  result  for  Ti  is  the 
same  as  the  result  of  Kubo  and  Tomita  theory. 

Nuclear  Spin-Lattice  Relaxation  for  Solid  Hvdrogen 
The  molecular  hydrogens  {H2,D2,HD,etc.)  form  the  simplest  molecular 
solids.  The  properties  of  solid  mixtures  of  ortho  (angular  momentum  J  =  1, 
nuclear  spin  /  =  1)  and  para  ( J  =  0,  /  =  0)  hydrogen  molecules  have  been 
extensively  studied  both  theoretically  and  experimentally  in  the  past  decade 
(53,54,55,56).  A  popular  method  of  experimentally  probing  this  system  has 
been  through  nuclear  magnetic  relaxation  studies  (57,58,59,60).  The  relax- 
ation is  determined  by  the  orientational  fluctuations  of  the  molecules  which 
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is  in  turn  determined  by  the  EQQ  interaction  between  the  ortho  molecules. 
This  relaxation  rate,  which  is  a  consequence  of  the  intramolecular  nuclear  spin 
interactions,  is  given  by  (56,G1) 

jr  =  jn^c^jlM  +  d^lljiM  +  y  j|(2u;o)]}  (2.68) 

Where  c  denotes  the  constant  of  spin-rotational  coupling  and  d  that  of  the 
intramolecular  dipolar  coupling,  with  respective  values  of  113.9  and  57.7  khz. 
The  spectral  density  functions  J{m(jj)  are  taken  at  wq  and  2ujq  where  loq  is  the 
Larmor  frequency(51,p278). 

We  consider  two  regimes. 
The  High  Concentration  Regime 

A.  B.  Harris  (56)  calculated  the  spectral  functions  for  the  correlation  func- 
tions for  both  infinite  and  finite  temperatures.  He  used  a  high  temperature 
expansion  method  to  calculate  the  second  moment  and  obtained  good  agree- 
ment with  the  high  concentration  experiments  of  Amstutiz  et  al.  (62),  with 
regard  to  both  the  temperature  and  concentration  dependence  of  the  relax- 
ation time. 

Myles  and  Ebner  (63)  used  a  high  temperature  diagrammatic  technique, 
combined  with  a  simlple  method  of  impurity  averaging  over  the  distribution 
of  O-H2  molecules.  The  averaged  equations  were  then  solved  numerically  to 
obtain  the  spectral  functions  for  solid  H2  self-consistently  for  the  first  time. 
The  resulting  spectral  functions  were  used  to  compute  the  Ti  as  a  function 
of  the  ortho-molecule  concentration  and  this  was  shown  to  agree  well  with 
experiments  at  lOK  and  over  concentraiton  range  of  0.5  <  X  <  1.  They 
obtained  a  y/X  concentration  dependence  for  Ti,  which  was  in  agreement  with 
the  data  of  Amstutz  and  colleagues  (62)  for  X  >  0.5.  i 
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The  Low  Concentration  Reghiie 

The  low  concentration  regime  {X  <  0.5)  had  been  explored  by  Sung(64), 
A.B.  Harris(56),  Hama  et  al.  (65),Ebner  and  Sung  (66),  and  Ebner  and  Myles 
(67)  at  an  earlier  date.  Recently,  the  work  has  been  concentrated  on  X  <  0.5 
and  very  low  temperatures  (2'  <  400mA;),  which  details  we  will  discuss  in 
Chapter  3.  ''  '  : 

Sung  (64)  applied  the  high  temperature  statistical  theory,  developed  for 
paramagnetic  resonance  with  a  small  concentration  of  spins,  to  the  calculation 
of  the  angular  momemetum  correlation  functions  and  Harris  used  an  improved 
version  of  the  same  theory.    The  Ti  resulting  from  these  calculations  had  a 

5 

concentration  dependence  of  Xs,  which  was  in  agreement  with  the  data  of 
Weinhaus  and  Meyer  (61),  but  the  magnitude  of  Tj  obtained  in  this  way  was 
in  disagreenient  with  that  data. 

Hama  et  al.  (65)  developed  a  theory  which  was  capable  of  treating  the 
T  =  oo  correlation  functions  at  all  concentrations  and  which  gave  a  concen- 
tration dependence  and  magnitude  for  Tj  which  were  in  fair  agreement  with 
experiment  for  all  X  (61,62). 

Both  methods  had  the  defect  that  the  impurity  averaged  correlation  func- 
tions were  obtained  by  statistically  averaging  assumed  functional  forms  and  no 
attempt  was  made  to  determine  the  shape  of  the  spectral  function. 

The  first  attempt  in  the  small  X  region  to  calculate  the  high  temperature 
correlation  functions  self-consistently  and  thus  to  overcome  the  above  defect 
was  made  by  Ebner  and  Sung  (66).  They  used  the  Sung  and  Arnold  (68) 
method  of  impurity  averaging  the  Blume  and  Hubbard  (69)  correlation  func- 

5 

tion  theory  and  obtained  a  Ti  which  had  the  experimentally  observed  Xs 
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concentration  dependence  at  small  X.  Since  they  made  no  attempt  to  prop- 
erly account  for  the  anisotropy  of  the  intermolecular  interactions,  they  did  not 
obtain  quantitative  agreement  with  the  experimental  magnitude  of  Ti. 

Ebner  and  Myles  (67)  improved  the  calculation  of  Ebner  and  Sung  by 
properly  treating  the  anisotropy  of  the  electric  quadrupole-quadrupole  (EQQ) 
interaction,  which  is  the  dominant  orientationally  dependent  interaction  be- 
tween two  O  —  H2  molecules  in  solid  H2  and  which  therefore,  almost  totally 
determined  the  shape  of  the  angular  momentum  spectral  functions. 

Sung  and  Arnold's  method  of  impurity  averaging  the  infinite  temperature 
Blume  and  Hubbard  (69)  correlation  function  equations  was  employed,  but 
the  equations  were  obtained  using  the  full  EQQ  interaction  rather  than  an 
isotropic  approximation  to  it.  The  spin-lattice  relaxation  time  was  computed 
as  a  function  of  the  O  —  H2  concentration  using  a  formula  for  y'  derived  by 
applying  the  Blume  and  Hubbard  theory  (69)  to  the  nuclear  spin  correlation 
functions  in  this  system.  The  resulting  Tj  was  compared  to  the  data  of  Wein- 
haus  et  al.  (61)  at  a  temperature  of  T  =  lOK  and  agreement  was  generally 
good  with  regard  to  both  its  magnitude  and  concentration  dependence. 


CHAPTER  3  .-.    .     ' 

I        •     . 

NUCLEAR  SPIN  LATTICE  RELAXATION 

Formulation  of  Longitudinal  Relaxation  Time  Ti 
There  is  a  striking  resemblance  between  the  phase  diagram  for  the  magnetic 
alloys  such  as  CuMn,  AuFe  ...    and  orientationally  ordered  ortho-hydrogen- 
para- hydrogen  alloys  (Fig.  1-2). 

The  nuclear  spin-lattice  relaxation  of  ortho  molecules  at  low  temperatures 
is  determined  by  the  modulation  of  the  intramolecular  nuclear  dipole-dipole 
interactions  Hjyjy  and  the  spin-rotational  coupling  Hgji  by  the  fluctuations  of 
the  molecular  orientations(55,  56).  The  calculations  are  particularly  transpar- 
ent if  we  use  orthonormal  irreducible  tensorial  operators  02m  and  A^2m  for  the 
orientational  [J  =  1)  and  nuclear  spin  (/  =  1)  degrees  of  freedom,  respectively. 

The  O2M  ^^^  given  by 

(3  J2  -  2) 

''^^/^ 

02±1  =  T-{J±Jz  +  JzJ±) 

02±2  =  \  {J±f  (3.1) 

and  similar  expressions  hold  for  the  N2m  in  the  manifold  7  =  1. 

The  intramolecular  nuclear  dipole-dipole  interaction  Hdjj  and  the  spin- 
rotation  interaction  IIsR  can  be  written  in  the  above  notation  as 


HDD  =  hD^Nljt)02m{i) 
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I 

I 


23    ;  i 

and  { 

HsR  =  -hcY,{-irN\^(i)0^m{i)  ;  (3.2) 

i 
respectively.     D=173.1  khz  and  C=113.9  khz.     The  index  i  labels  the  ith 

molecule.    The  0;„,  and  A^;^  are  the  operator  equivalents  of  the  spherical 

harmonics  F/^  in  the  manifolds  J  =  1  and  /  =  1,  respectively. 

The  relaxation  rate  due  to //i>p  can  be  shown  to  be  j 

=  ;^^"([/.,>/i)Z?]  [er^^''>iDDe-^"'\lz])dt  .       (3.3) 

where  Ho  —  hojoh-   It  is  the  Hamiltonian  responsible  for  the  molecular  dy- 
namics. Using  the  commutators  [Iz,N2m\  ~  i^^2m  we  find 


'^IDD  =  \d'^     E     rn^J2m{mu^o)  i  (3-4) 

|m|  =  l,2  ■ 

where  the  spectral  density  at  the  Larmor  frequency  .  i 

/oo  .  . 

(O2m(0oL(0))T«~"'°*^^  ;  '.      (3-5) 

-oo 

The  expectation  value  <  ...  >t  must  be  calculated  with  respect  to  the 

fixed  Z  axis  given  by  the  direction  of  the  external  magnetic  field.    This  is 

the  quantization  axis  for  the  nuclear  Zeeman  Hamiltonian,  which  is  perturbed 

by  the  weaker  H[)d  and  Hg^i  terms.    The  orientational  order  parameters, 

however,  are  evaluated  with  respect  to  the  local  molecular  symmetry  axes.  We 

must  therefore  consider  the  rotations 

02m  =  E^S(x)02/*  ;       -         (3.6) 

where  the  tirn/x  are  the  rotation  matrix  elements  for  polar  angles  x  =  (q:>/^) 
defining  the  orientation  of  OZ  in  the  local  (x,y,z)  reference  frame. 
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Figure  1.  Phase  Diagram  of  Ortho-Para  H2  Mixtures 
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(A)  LONG  RANGE 
ORDER 


(B)  GLASS 


Figure  3.  Picture  of  Long  Range 
Order  And  Quadrupolar  Glass 


27 


We  assume  the  simplest  possible  case 


{02m{t)oljO))  ={O2m{0)olj0))g2m{t) 


(3.7) 


In  the  following,  we  consider  only  this  case  and  further  assume  that  the  relax- 
ation is  dominated  by  the  fluctuations  of  the  /Lt  =  0  component.  We  find 


^ll!z?  =  r2(^-^-^^)^^    E    "^^\dmoH\^92o{mujQ)  (3.8) 

m=l,2  , 


where  the  g2o{''^^o)  are  the  Fourier  transforms  of  the  reduced  correlation  func- 
tions §20(0  ^^^  ^  =(3'^z  -  2)2-.  The  prefactor  (2  -  a  -  a'^)  is  the  mean  square 
deviation  of  the  operator  O20  evaluated  in  the  local  symmetry  axis  frame. 
The  contribution  from  the  spin-rotational  interaction  Hgji  is 

T[sR  =  l'^''jnM  =  lc^{2  +  o)\d,o{a)\^g,oM  ^    (3.9) 
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the  total  rate  T^  ^  =  T^^J^  +  T~gj^ 


Temperature  Dependence  of  Ti 


Minimum  Values  of  Ti 


From  the  previous  result  (3.8  and  3.9) 


where 


^1  ^IDD  ^ISR 


—  {o,e)  =  ^{2-a-a^)D^    J2    m^\dmOm^g2o{muo) 


Two 


m=l,2 


If  we  take  a  powder  sample  average,  the  Tigji  is  small. 


^W=.l,2-.-.V^ 


1  A       ^         4  , 

^ff2o(wo)  +  gff2o(2wo) 


(3.10) 
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To  a  good  approximation  we  have 

^(a)«^(2-a-a>2 


(3.11) 


For  the  simplest  case  we  can  restrict  a  to  negative  values  with  appropriate 
definitions  of  principal  axes  (22  and  Chapter  4) 


I   _/°,(2-a-a2)da      j_ 


1 


7-2  <^^ 

1 


12 


g&2o(woj  +  -ff2o(2a;o) 


=  ^^'(^2o(wo)+4y2o(2u;o)) 


(3.12) 


Assume  fif2o('*'o)  and  ff2o(2wo)  can  be  taken  to  be  in  the  Lorentzian  form: 


920  (i^o)  = 


Tc 


ff2o(2wo) 


l  +  u;2^2 
1  +  4a;2^2 


(3.13) 
(3.14) 


when  UQTc  «  0.6156,  Ti  =  T^^i^ 

The  eqs.  (3.12)  —  (3.14)  result  in  the  following 

1) 

Wo  =  27r  X  100  X  10® 

^Imin  —  13.4172  msec 


2) 


Wo  =  27r  X  25  X  10** 

Timin  =  3.3542  msec 

Tuxp  =  1.03  msec  -  Fig  A  (ref.25) 


3) 


Wo  =  27r  X  9  X  10** 

^Imin  ~  1-20574  msec 

Tiexp  =  2.25  msec  -  Fig.5  (ref.70) 
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Figure  5.  Experimental  Curve  of  Ti 
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Temperature  Dependence  of  Ti 

It  has  been  observed  that  for  the  solid,  a  Gaussian  Free  Induction  Decay 
is  a  good  approximation  for  small  t  and  so  we  should  also  consider  a  Gaussian 
form  for  g-zoirnujQ)  for  high  frequencies.  From  eqs.  (3.7)  and  (3.8): 

{O2m{0)olj0))g{t)  =(O2m(0)O^^(0))e-S         ;  (3.15) 

-00 
Eq.  (3.12)  becomes 


—- =  ^^—(ae       4     -\-4ae   '^o"  i 

Ti  36 


(ae      4     +  4ae~'^oa  )  .  (3.17) 


^  36  1 

Ti  =  —2-7= Y^ (3.18) 

Ti  passes  through  a  minimum  when  ^  =  0,  i.e  for  a  =  5.082  x  10~^,  and 
the  minimum  value  is  Tj^^^  =  1.1384msec,  at  |^  =  25  MHz.  The  result  of 
the  theoretical  value  T^y^j^  =  1.1384m5ec  is  in  very  good  agreement  with  the 
experimental  value  Ti^i^  =  1.03  msec  (25). 

Now  the  question  is  how  Tj  varies  at  fixed  uq  over  a  wide  range  of  tem- 
peratures which  causes  Tc  to  vary.  Similar  to  the  discussion  of  the  dynamics  of 
spin  glasses  (71),  we  would  like  to  try  an  Activation  Law  of  the  Vogel-Fulcher 
type: 

a  =  aot^^^'^l  (3.19) 

where 

A  is  the  activation  energy  in  temperature  unit 

oq  is  a  time  factor  which  is  the  value  of  a  as  T"  ^^  00,  and 

To  is  some  characteristic  temperature  (transition  temperature),  such  that 

as  r  — >  Tq,  long  relaxation  times  become  important. 
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Figure  6.  Ti-T  Curve  (solid-TiexpA^iied-Titheo) 
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By  using  two  experimental  values  (Ti  =  2.550msec  at  T  =  0.450i^  and 
Ti  =  2.514msec  at  T  =  0.380K  for  X  =  38%)  we  find  the  following  formula: 

_o  0.1880 

a  =  5.4128  X  10   °  X  e^-o.oois  (3.20) 

The  Fig.6  shows  curves  of  results  of  calculation  and  experiment.  The  solid 
line  is  experimental  curve(25)  and  the  dashed  line  is  theoretical  curve. 

Spectral  Inhomogeneitv  of  Ti 
Results  of  Calculations 

The  dependence  of  Ti  on  Au  was  evaluated  by  considering  the  line  shape 
to  be  a  sum  of  Pake  doublets.  Each  Pake  doublet  consists  of  a  positive  branch 
given  by  Av  =  +^DP2{a)a  and  a  negative  branch  given  by  Av  =  —^DP2{a)a 
(Fig. 7).  As  previously  demonstrated  we  restrict  a  to  negative  values.  That  is 

Au  =  —D (3cos  a  —  1)  [positive  branch) 

Li      £i 

Au  —  D—-{Zcos  a  —  1)  {negative branch) 

Ai     At 

In  order  to  test  the  theory  against  the  experimental  data  we  consider  only 
the  low  temperature  limit  for  which  the  fluctuations  are  slow  compared  to  the 
Larmor  frequency  wq.  In  this  case  the  spectral  density  functions  are  given  by 
S'('^o)  =  4ff(2wo)  =  -5-  (low  temperature  limit  of  Lorentz  form),  where  r~^  is 
the  characteristic  fluctuation  rate,  which  was  taken  to  be  a  unique  value  for 
simplifying  calculations. 

For  the  dipolar  contribution  at  fixed  o  and  Au  we  obtain 


rr-\      _ 

A(2- 

-  0  - 

-  a^)D^ 

1  - 

(2Aiy\2 

w^r 

(3.21) 
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Figure  7.  Allowed  Range  of  Values  of  o  for 
A  Given  Frequency  Ai/ 
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and  for  the  spin-rotational  contribution: 

-1C2(2-H)(1±^) 


ISR 


u;2r 


(3.22) 


The  frequency  dependence  should  be  obtained  by  summing  over  all  allowed 
a  for  a  given  frequency  Ai/.  The  calculations  were  straightforward  but  quite 
tedious.  Table  I  and  Table  II  sPiow  the  numerical  results. 


Table    I    Spectral    Dependence   of   Dipole-Dipole    Relaxation    Rate    T.J; 


IDD 


Frequency  range  Ai/ 

0<Au<^D 

WD''-12uD+4u^    1         2i^-'       ,_2i/    ,        2i.'^       .^u 

3£>2                1    D{D-2u)  '°e  D    1    DiD-i.)  '°S  D 

\D  <Au<D 

BD'-luD+iu^    I        2u'       .    ^  u 

3Z?2              1    D(D-u)  '°S  D 

Table  II  Spectral   Dependence  of  Spin-Rotational  Relaxation   Rate   T,  ^ 


ISR 


Frequency  range  Au 

0<Au<Id 

D            D-2v  '«g  D    1    D-v  'Og  V 

\D  <Au  <D 

D-Vv    1      2i/     i„„  V 
D     1   £>-<.Jogn 

Discussion 

Fig. 8  shows  the  results  of  the  calculations  for  the  dipolar  contribution  TrJj^ 
and  the  spin-rotational  contribution  T^gji  as  function  of  frequency.  All  curves 
in  this  figure  have  been  normalized  to  unity  at  Au  =  0  in  order  to  facilitate 

the  comparison  with  the  experimental  data.  The  net  relaxation  rate  T7i^  ,  <  = 

l(calc) 

^IDD  "^  '^ISR'  '^^^  curves  show  discontinuities  in  slope  at  |Ai/|  =  y,  but  this 
was  not  seen  experimentally.  If  we  take  the  finite  cross-relaxation  rate  TTi}  into 


':■-:'«* 
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account,  which  will  bring  the  very  slowly  relaxing  isochromats  at  |Aj/|  =  2D 
into  communication  with  the  rapidly  relaxing  components, 

2^r'  -  ^iMc)  +  ^1-2'  (3.23) 

■.     •    -     ■■  .1 

when  the  individual  spin-lattice  relaxation  times  are  much  longer  than  the 
cross-relaxation  times  the  spins  come  to  a  common  spin  temperature  via  the 
cross-relaxation  mechanism  before  relaxing  to  the  lattice  via  the  rapidly  relax- 
ing components.  On  the  other  hand,  when  the  direct  spin-lattice  relaxation  is 
fast  and  T^  <  T12,  the  magnetization  of  each  molecule  relaxes  directly  to  the 
lattice  and  the  spins  do  not  achieve  a  common  spin  temperature.  In  this  latter 
case  the  nuclear  magnetization  will  be  spatially  inhomogeneous. 

Yu  et  al.  (14)  defined  Ty2  by  the  probability  of  a  spin  flip-flop  transition 
via  the  intermolecular  dipole-dipole  interaction  for  two  isochromats  ui  and  1^2 
given  by 


T^^  =  (r/'*>-'^'°P)-lexp 


7r^(/^l  -  U2)'' 
Jointer 


^Tf'P-f^<'P)--^exp 


jifintra  _  /  ii^intra\2 

2MY^^ 


(3.24) 


where  the  exponential  factor  is  the  overlap  given  by  Abragam  (51)  for  two 
lines  centred  at  fj  and  U2  and  their  individual  widths  are  determined  by  the 
intermolecular  dipole-dipole  interaction.  M^^^^^,  j^tntra  ^^^^  M^^^^''  are  the 
moments  resulting  from  the  intra-  and  inter-molecular  dipole-dipole  interac- 
tions, respectively.  Yu  et  al.  (14)  observed  a  much  weaker  dependence  than 
the  exponential  variation  with  M^^*^^  given  by  the  overlap  factor.  Due  to 
the  discrepancy  between  theory  and  experiment  we  chose  the  empirical  values 
reported  by  Yu  et  al.  (14)  as  the  most  reliable  estimate  of  T12.    . 


37 


500     - 


100     - 


T,(Z/J 
T,(o) 


10     - 


Figure  8.  Frequency  Dependence  of  Ti 
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The  most  complete  studies  of  the  spectral  inhomogeneity  of  Ti  in  the 
quadrupolar  glass  phase  has  been  carried  out  for  an  ortho  concentration  X  = 
0.45  at  T  =  0.15/<'(14),  and  the  data  of  ref.l4  would  place  the  cross-relaxation 
time  Ti2  in  the  range  7.5  —  10.5  msec. 

While  the  cross-relaxation  is  faster  than  the  direct  relaxation  to  the  lattice, 

the  values  of  2^2  reduce  the  spectral  inhomogeneity  of  the  relaxation  of  the 

NMR  line  shape.  In  this  case  we  expect  to  observe  TT,  ,  >  =  T~,    ,.  +  TTo    for 
^  ^  l(o6s)  l(cal)  1'2 

the  relaxation  of  the  magnetization  of  a  given  isochromat.  This  is  shown  in 
Fig. 8.  In  view  of  the  uncertainties  in  the  cross-relaxation  and  the  simplifying 
assumptions  that  have  been  made,  the  overall  agreement  with  the  experimental 
results  is  good.  The  correct  overall  behavior  is  predicted  as  well  as  the  subtle 
change  in  spectral  dependence  at  the  half  width  points  which  has  already  been 
seen  in  the  experimental  data  of  Yu  et  al.  (14).  ' 

Non-Exponential  Relaxation  of  Nuclear  Magnetization 
Although  it  is  agreed  that  (i)  the  low  temperature  NMR  lineshapes  indi- 
cate a  random  distribution  of  molecular  orientations  (for  both  the  local  axes 
and  the  alignment  a  =(3J^  ~  2)),  and  (ii)  that  there  is  apparently  no  abrupt 
transition  in  the  thermodynamic  sense;  there  has  been  disagreement  (1,3)  over 
the  interpretation  of  the  behavior  of  the  molecular  orientational  fluctuations 
on  cooling  from  the  high  temperature  (free  rotator)  phase  to  very  low  tem- 
peratures (r  <  O.IK).  Is  it  a  very  rapid,  but  smooth  variation  (25,27,70) 
with  temperature,  corresponding  to  a  collective  freezing  of  the  orientational 
fluctuations  or  a  slow,  smooth  dependence  (9,14,28)  with  no  evidence  of  any 
strong  collective  behavior?    In  order  to  resolve  this  problem  it  is  important 
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to  understand  two  unusual  properties-tlie  spectral  inhomogeneity  of  the  re- 
laxation rate  T^  across  the  NMR  absorption  spectrum  (14,29,30,72)  and  the 
non-exponential  decay  of  the  magnetization  of  a  given  isochroinat  (14).  The 
spectral  inhomogeneity  has  been  discussed  previously  and  the  purpose  of  this 
part  is  to  show  that  the  strong  departure  from  exponential  decay  for  the  mag- 
netization can  be  understood  provided  that  the  broad  distribution  of  local 
order  parameters  is  correctly  accounted  for. 
Variation  of  Relaxation  Rates  within  Given  Isochromats 

As  previously  proved  (3.21,3.22)  I 


2 


k{2-a-a^)D^    1  -  (^) 


24 


,-1  

IDD  -  ,  ,2, 

2At/> 


wgr 


7^-1    _ 

-I  1  O  D    — 


l(2-|a|)C2(l±^) 


18 


The  frequency  of  a  particular  component  of  the  NMR  absorption  line  is 
given  by  • 

Au  =  ±-DP2{a)a  (3.25) 

and  this  can  be  satisfied  by  very  different  values  of  a  and  a;  e.g.  Aw  —  ^D 
occurs  for  a  =  -1,^2  =  |;  <^  =  -|»^2  =  1;<^  =  -§)  ^2  =  f---;  the  only 
constraints  being  that  a  and  P2  lie  within  their  limits;  — 2  <  a  <  0,  and 
—  ^  <  ^2  ^  1-  Obviously,  different  pairs  of  a  and  P2  (for  fixed  Au)  can  result 
in  very  different  values  of  T^pp  and  TisR-  Molecules  which  contribute  to  the 
same  isochromat  of  the  NMR  line  but  which  have  different  values  of  a  and  P2 
will  therefore  relax  at  different  rates.  This  is  illustrated  in  table  III,IV  and 
V,  which  give  the  variation  of  Ti  for  Au  =  \D,  \D  and  0,  respectively.  The 
relative  contribution  of  these  rates  can  be  determined  from  the  probabilities 
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n(a),n(P2)  of  finding  a  and  P2.  At  low  temperatures,  the  analysis  of  the  line- 
shape  indicates  that  a  good  approximation  for  IT  (a)  is  a  triangular  distribution 
n(a)  a  o  (12,23).  For  P2  we  assume  a  powder  distribution  of  local  axes  (i.e. 


very  glassy)  which  requires  that  n(P2)   c< 


)/^" 


The  calculated  rates 


Tj"  (cr,  P2)  have  a  relative  weight  P  =  Ti[a)  x  n(P2)  for  the  relaxation  of  the 
component  t\u  =  '^^D\o\P2. 
Magnetization  ratio 

■  m__ZP^  .326) 

M(0)  "       Y.P  ^  ^     ^ 

The  weighted  relaxations  using  the  indicated  probabilities  are  given  in  Fig. 9 
and  Fig. 10.  The  experimental  results  reported  by  Yu  et  al.  (14)  for  different 
/\u  are  indicated  by  the  symbols. 

Table  III  Variation  of  Relaxation  Rates  within  A  Given  Isochromat  (Ai^  =  ^D) 


TatesI^ 


parame. 

IF 
ilj- 


parame. 

kl 

_15 


prob. 

n(P2) 


prob. 

n(a) 


rates 


rates 


1.022 


1.067 


0.058 


0.234 


0.267 


-if 


1.044 


1.143 


0.107 


0.431 


0.492 


-1^ 


1.069 


1.231 


0.144 


0.583 


0.666 


J3- 


1.095 


1.333 


0.167 


0.681 


0.777 


1.124 


1.455 


0.170 


0.706 


0.804 


T 
-1- 


-/- 


1.155 


1.600 


0.150 


0.634 


0.720 


1.188 


1.778 


0.097 


0.422 


0.478 


1.225 


2.000 


0.0 


0.0 


0.0 


("J Rates  given  in  units  of  -T  2 


The  Determination  of  The  Molecular  Correlation  Time  t  and  T^  at  Ai/  =  0 

The  only  unknown  parameter  for  eqs.    (3.21)  and  (3.22)  is  the  molecular 
correlation  time  r  and  the  best  fit  represented  by  the  solid  lines  and  the  broken 
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Table  IV  Variation  of  Relaxation  Rates  within  A  Given  Isochromat  (Ai/  =  |Z?) 


parame. 
IS" 


3arame. 


f 


"X" 


-f 


-13- 


_LL 


prob. 
W2) 


1.022 


1.044 


1.0G9 


1.119 


1.155 


1.188 


1.225 


1.265 


1.309 


1.359 


1.414 


prob. 


0.533 


0.571 


0.615 


0.667 


0.727 


0.800 


0.889 


1.000 
1.143 


1.333 


1.600 


2.000 


rates 

rp-1 


0.092 


0.179 


0.260 


0.333 


0.398 


0.450 


0.486 


0.500 


0.482 


0.417 


0.275 


0.0 


rates 

rp-l 

JLiniL 


0.272 


0.526 


0.760 


0.972 


1.159 


1.316 


1.435 


1.500 


1.485 


1.337 


0.938 


0.0 


rates 


R 


^-1 


0.325 


0.629 


0.910 


1.165 


1.389 


1.576 


1.715 


1.789 


1.763 


1.577 


1.097 


0.0 


v**' Rates  given  in  units  of  „P ^ 

24u}^T 


line  at  Az/  =  0  in  Fig.9,  is  obtained  t  =  1.51  x  lO^'^S'.  For 


T  =  1.51  X  10~'''5 


1 


10"^ 


^   Az/=0=  1-49597  X 

Ti  '^"^    "  10.0446 

Ti  |ai/=o=  6.7l44(msec) 


(3.27) 


Ti  at  Ai/  =  0  is  6.71  ms  which  is  in  excellent  agreement  with  the  experimental 
value  of  6.6  ±  0.5ms  (14,70). 

Since  D=173.1  khz,  in  Fig.9  the  Au  for  the  calculated  ^14  are  87,43.5  and 
0  khz  which  are  below  the  values  chosen  by  the  Duke  group  (14).  The  same 
theroy  to  calculate  j^  for  Ai/  =  98  and  58  khz  and  the  results  are  depicted 
in  Fig. 10.  The  overall  agreement  is  very  satisfactory. 

Comparison  of  the  calculated  decays  M^j^{t)  with  the  experimental  results 
shows  that  not  only  is  the  correct  overall  deviation  from  exponential  decay 
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predicted,  but  that  thei'e  is  also  a  significant  long  tail  to  the  decay  which 
ought  to  be  tested  for  experimentally.  This  long  time  behavior  is  unique  to  the 
glassy  regime  of  the  hydrogen  mixtures. 


Table  V  Variation  of  Relaxation  Rate's  within  A  Given  Isochromat  {Au  =  0) 


paranie. 

kl 


-15. 


-f 


^ 


prob. 
11(a) 


0. 


0.571 


0.615 


0.667 


0.727 


rates 


1.9/? 


1.467 


1.429 


1.385 


1.333 


1.273 


rates 

rp-l 


2.249 


2.245 


2.237 


2.222 


2.198 


rates 


R" 


.-1 


3.096 


3.070 


3.036 


2.992 


2.933 


0.800 


1.200 


2.160 


2.853 


0.889 


1.111 


2.099 


2.740 


1 


1.000 


1.000 


2.000 


2.577 


1.143 


0.857 


1.837 


2.331 


1.333 


0.667 


1.556 


1.940 


1.600 


2.000 


0.400 


0.0 


1.040 


0.0 


1.271 


0.0 


We  consider  P2(«)  =  0  for  the  range  of  \a\  considered  in  Table  IV  in  order 
to  facilitate  the  comparison. 

''^^ Rates  given  in  units  of  ;r^ 


24u;2r 
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Figure  9.  Time  Dependence  of  M(t)  for  Different  Ai/ 
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Figure  10.  Time  Dependence  of  M(t)  for 
Au  =  98  and  Au  =  58  khz 


CHAPTER  4 
ORIENTATIONAL  ORDER  PARAMETERS 

Density  Matrix  Formalism 
A  particle  (e.g.,  an  atom,  molecule  or  nucleus)  isolated  in  space  and  with 
nonzero  angular  momentum  in  ito  rest  frame  has  a  manifold  of  states  with  equal 
energy.  The  problem  is  how  to  specify  the  orientational  degrees  of  freedom  of 
an  individual  molecule.  We  describe  the  degrees  of  freedom  of  the  quantum 
rotors  with  angular  momentum  J  =  1  in  terms  of  single  particle  3x3  density 
matrices  p^  (for  each  site  i).  The  p^  are  completely  described  by 

(1)  the  molecular  dipole  mements  {Jx)^,  {Jy)i,{Jz)i  and 

1 

(2)  the  quadrupole  moments  (J^)^,  {JxJy)i,  {JyJz)^....  In  the  absence  of 
interactions  which  break  time  reversal  symmetry,  the  dipole  moments  {Jx), 
(Jy),  and  (Jz)  vanish  and  we  need  only  consider  5  independent  variables. 

Instead  of  Cartesian  components,  it  is  more  convenient  to  use  a  set  of 
irreducible  tensorial  operatiors  Ulm  with  0  <  L  <  2J  for  general  J  and  the 
associated  multipole  moments 

ti^  =  TripUi^)  (4.1) 

For  simplicity  the  site  index  has  been  dropped.   The  expansion  of  the  single 

particle  density  operator  in  terms  of  the  multipole  moments  is  given  by 

2J       L 

f=J^7TT)T.     E    i'^L  +  l)tlj^ULM  (4.2) 

^  '  L=--OM=-L 

There  are  three  conditions  imposed  on  p:  Hermiticity  and  both  weak  and 
strong  positivity  conditions(34). 

45 
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/)  is  a  HcrmiLlan  operator  and 


\Mj* 


tLM=i-irtl-M  Oi-3) 


Uoo  is  a  unit  matrix  operator  and  too  =  1  =  Trp 
The  weak  positivity  conditions  are  given  by 

27117  ^ '''•''''*  ^ '  ;  <'-^' 

When  one  eigenvalue  Js  equal  to  1,  the  others  being  null;  then  the  matrix  p 
describes  a  pure  state  and  Tr{p'^)  =  1.  The  minimum  of  Tr{p'^)  is  reached 
when  all  the  eigenvalues  are  equal  to  (2J  +  1)~^;  then  the  matrix  describes  a 
completely  unpolarized  slate  and  rr(/9'^)  =  (2 J  +  1)~^. 

,  For  density  matrices  of  spin— ^  particles,  condition  (4.4)  is  the  only  con- 
dition imposed  by  the  positivity  condition.  But  for  J  >  \,  further  conditions 
are  imposed  on  the  density  matrix  and  on  the  multipole  parameters  by  the 
positivity  property.  ,     . ,  ,         ,  . 

The  eigenvalues  of  /?  must  be  positive  definite  because  they  represent  the 
probabilities  of  realizing  some  given  state  and  this  leads  to  the  strong  positivity 
conditions 

0  <  An  <  1  '  (4.5) 

where  An  is  the  nth  eigenvalue.  These  conditions  place  the  strongest  limitations 
on  the  allowed  values  for  the  multipole  moments  and  thus  on  the  allowed  values 
of  the  local  order  parameters  for  ortho-^2  molecules  in  the  solid  mixtures. 

It  is  useful  to  construct  orthonormal  matrix  representations  of  the  irre- 
ducible operators  IT^jj/  in  the  representation  ( J'^,  J^).  For  J=l  these  are  given 
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by  the  following  3x3,  iuat<ix  operators  with  rows  (and  columns)  labelled  by 
the  eigenvalues  1,0,-1  of  Jz. 

1  1    /I     0      0 

.    !■  '    .    .V?    . ,    v2  Vo    0    -1 

1  1    /O    1    0- 

IIll  =  — ^J+  = ^     0     0     1 

v2  V2  \o     0     0, 

n20  =  ^(34-^')  =  ^(0     -2     OJ  (4.6) 

1  1    /o     1      0 

■   ^21  =  --(./^J+  M-  J+J;^)  =  -^     0     0-1 

-^  v2  \o     0       0 

1  /o    0    i\ 

.         .        ,n22  =  ;^4=      0     0     0  (4.7) 

^         Vo   0   0/ 
ni,M  =  (-i)''nt    ^  (4.8) 


and 


Tr{TlLM^\>M')  =  ^LL'^MM'  (4.9) 


The  quadrupole  operators  n2M  transform  analogously  to  the  spherical  har- 
monics Y2m{oi,l3)  with  respect  to  rotations  of  the  coordinate  axes. 

The  reference  axes  have  remained  arbitrary  in  the  discussion  and  we  are 
therefore  free  to  choose  local  references  axes  that  correspond  to  the  local  sym- 
metry for  each  molecule.  The  natural  choice  for  the  z-axis  is  along  the  net  com- 
ponent of  the  angular  momentum  at  a  given  site,  i.e.  such  that  (Jx)  ={Jx)  =  0. 

The  general  form  for  /Jj=i  can  be  identified  by  the  mean  values  of  its 
magnetic  dipole  and  electric  quadrupole  rnements  in  the  above  notation: 

^=  on3  +  yi^'^n„  + Vg,„n2^  (4.10) 

n  m 

where  /i„  =z(ri[^)  and  Qm  =(^1^) 
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The  expression  (4.10)  is  identical  with  that  was  derived  by  ref.35  and  ref.73. 
We  still  can  choose  x  and  y  axes  such  that  Q2  is  real,  i.e. 


1^ 


Q2  =  ^{4-J^) 


(4.11) 


and 


{JxJy  +  JyJx)  =  0  (4.12) 

Q2  measures  the  departure  from  axial  symmetry  about  the  z  axis  and  is  some- 
times called  the  eccentricity  (29,  30).  It  can  be  shown  that  with  this  choice  of 
local  reference  frame  Qi  and  Q_i  also  vanish  and  the  density  matrix  may  be 


written  as 


1  1         /I     0      0 

p  =  -Hs  +  —=^lQ    000 

^  v2       Vo     0 


/^fe^o 


where 


0      0       + 

0    -1/ 

0 

[     Q2 

/^o=^W 

C?o  =  >3 

Jl  -  2)) 

v/6 


Qo 


Q2    \ 

0 


(4.13) 


(4.14) 

Qq  is  the  alignment  (29,  30)  along  the  z-axis.    Sometimes  it  is  convenient  to 
define  the  alignment  and  the  eccentricity  as 


o 


il--,J^)- 


Iqo 


(4.15) 

7?  =(4  -  j2)  =  2Q2  (4.16) 

Both  a'  and  7/  have  maximum  amplitude  of  unity.  In  terms  of  these  parameters 


p  becomes 


1     0      0 


p  =  -Hg  +  -//     0     0 


+ 


0     0-1 


1^' 

-3^ 

0 

in\ 

0 

i"' 

0 

b 

0 

-1-7 

(4.17) 
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/i*+i?*=j/ 


Figure  11.  The  Allowed  Values  of/i  ={Jz), 
a  And  77  for  Spin-1  Particle 


m 


50 

.     •  i 

'     '  '  i 

where  /x  ={Jz). 

The  three  eigenvalues  of  p  are 

12, 
^^  =  3  +  3^^ 

A2,(3)  =  ^-^^'±^\^^  +  ^  '  ('-'^^ 

The  strong  positivity  conditions,  An  >  0,  are  therefore  seen  to  restrict  the 
allowed  values  of  the  local  order  parameters  a' ,r]  and  /x  to  the  interior  of  a 
cone  (Fig.  11)  in  the  3D  parameter  space.  The  vertex  of  the  cone  is  located  at 
a'  =  1,  fj,  =  t]  =  0,  corresponding  to  the  pure  state  \  ip)  =\  Jz  =  0) ,  and  the 
base  of  the  cone  is  defined  by  a'  =  -^  and  /x^  +  r/^  =  1  which  corresponds 
to  the  pure  state  |  tp)  =  cos'^  \  Jz  =  1)  +  sin^  |  Jz  =  -1)  with  (Jz)  =  cos2^ 
and  ( J^)  =  sin2^  (The  polar  angle  7  generates  the  points  on  the  circle  of  the 
cone's  baseplate). 

The  positivity  domain  shown  in  Fig. 11  is  the  same  as  those  obtained  by 
Minnaert  (34)  using  the  Ebhard-Good  theorem  and  similar  to  those  given  by 
W.Lakin  (32)  in  his  analysis  of  the  states  of  polarization  of  the  deuteron. 
Having  established  the  physical  considerations  which  determine  the  limited 
range  of  allowed  values  for  the  order  parameters,  we  now  turn  to  the  special 
case  of  solid  hydrogen. 

Applications  to  Solid  Hydrogen 
In  the  absence  of  interactions  which  break  time  reversal  symmetry,  the 
expectation  value  (Ja)  must  vanish  for  all  a  in  solid  hydrogen.  This  is  the 
so-called  "quenching"  of  the  orbital  angular  momentum  (74).  The  reason  for 
this  is  that  in  the  solid  the  electronic  distribution  of  a  given  molecule  may 
(to  a  first  approximation)  be  regarded  as  being  in  an  inhomogeneous  electric 
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field  which  represents  the  effect  of  the  other  molecules.  This  inhomogeneity 
removes  the  spatial  degeneracy  of  the  molecular  wave  function  which  must  be 
real  and  the  expectation  value  of  the  orbital  angular  momentum  {ih-ir)  must 
accordingly  vanish.  The  separation  of  the  rotational  energy  levels  is  given  by 
Ej  =  BJ{J  + 1)  with  B=85.37K,  and  in  the  solid  at  low  temperatures  only  the 
lowest  J  values,  J  =  0  for  para-/f2  '^'^^  J  =  1  for  OTtho-H2,  need  be  considered. 

At  low  temperatures  the  anisotropic  forces  between  the  ortho  molecules  lift 
the  rotational  degeneracy  and  the  molecules  align  themselves  with  respect  to 
one  another  to  minimize  their  interaction  energy.  For  high  ortho  concentrations 
one  observes  a  periodic  alignment  in  a  Pa^  configuration  with  four  interpene- 
trating simple  cubic  sub-lattices,  the  molecules  being  aligned  parallel  to  a  given 
body  diagonal  in  each  sub-lattice  and  the  order  parameters,  <7^  =(l  —  |j|a)» 
are  the  same  at  each  site  a^  =  1. 

The  long  range  periodic  order  for  the  molecular  alignments  is  lost  below 
a  critical  concentration  of  approximately  55%  and  the  NMR  studies  indicate 
that  there  is  only  short  range  orientational  ordering  with  a  broad  distribution  of 
local  symmetry  axes  and  local  order  parameters  throughout  the  sample.  This 
purely  local  ordering  has  been  referred  to  as  a  quadrupolar  glass  in  analogy 
with  the  spin  glasses,  but  unlike  the  dipolar  spin  glasses,  there  is  no  well-defined 
transition  from  the  disordered  state  to  the  glass  regime.  In  order  to  describe 
the  ortho  molecules  in  the  glass  regime,  where  there  is  a  large  number  of  sites 
with  various  values  of  a',  a  density  matrix  formalism  must  be  used. 

The  quenching  of  the  angular  momentum  in  solid  H2  has  two  consequences 
for  the  limits  on  the  allowed  values  for  the  order  paramters:         '  ' 
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(1)  From  the  previous  discussion,  the  allowed  values  of  o'  and  rj  lie  within 
a  triangle  bounded  by  the  three  lines  (Fig.  12  l\ABC).  \. 

,-12,  •    r    • 

-  +  -cr'  >  0 
.....  3       3      - 

,     ^--a'±-r/>0         ,      ,  .   .      (4.19) 

which  represent  the  strong  positivity  conditions  for  the  eigenvalues  of  p  when 
the  angular  moinentuui  is  quenched.  , 

(2)  Since  {J^)  —  0  for  all  a,  we' are  free  to  choose  the  z-axis  which  was 
previously  fixed  by  the  net  component  of  the  angular  momentum.  The  natural 
choice  for  the  local  reference  axes  is  now  the  set  of  principal  axes  for  the 
quadrupolar  tensor  '  , 

.       .  ,  ..        Qafi={\[JaJl3  +  JpJa)-\j'^B^^)  '  (4.20) 

The  choice  of  principal  axes  is  not  unique,  however,  because  after  finding  one 
set  we  can  always  find  another  five  by  relabeling  axes.  We  can  prove  that  not 
all  of  the  points  in  the  "allowed"  portion  of  the  [o' ,r})  plane  are  inequivalent 
(every  point  represents  one  state),  and  we  only  need  to  consider  the  hatched 
region  (triangle  CFE)  in  Fig. 12.  The  states  of  all  other  areas  of  the  triangle 
of  allowed  values  can  be  obtained  from  the  ACFE  by  suitable  rotations  (i.e. 
relabel  the  axes). 

The  important  point  is  that  the  pure  states 

^C=\Jz  =  0)  ,  (4.21) 

and 

^'A  ----  -kl  -^^  =  1)+  I  J^  =  -i))  H-22) 

V2 
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are  not  inequivalent  (The  labels  A  through  F  refer  to  the  special  points  on 
the  triangle  of  allowed  values  shown  in  Fig. 12),  and  the  corresponding  wave 
functions  are  listed  in  Table  VI.  The  state  V'A  can  be  obtained  from  xl>c  by  a 
rotation  of  the  axes  by  ^-  about  x  axis. 
The  rotation  operator 

•         ■    '  ••       RA~)  =  ^xp{-^~J.)  =  -iff  +  ^{f)  +  I    ■  (4.23) 

and 

I^x[^)^C  =  e'2^^  (4.24) 

The  rotation  Rx{-y)  leaves  the  state  F  invariant  and  maps  D  onto  the  point 
G  in  parameter  space. 

We  can  furthermore  show  that  the  rotation  Rx{^-)  maps  the  following 
triangular  regions  of  parameter  space  onto  one  another 

ACPF  ^  IsAGF 

ACEF  ->  AAEF 

and 

ABDF  -^  ABGF  '  (4.25) 

It  is  also  seen  by  considering  the  transformation 


/(a",  V)  =  4(y)p((r',;7)ii:.(y)  (4.26) 


using  the  matrix  representation 


«.(f )  = 


/     2         ^  2    \ 

v/2      "      v^ 

V      2       ^        2    J 


(4.27) 


54 

Table  VI  Special  Points  m  Oiieutalional  Order  Parameter  Space,  (see  Fig.  12) 


Label 

A 

paraine. 
(-ll) 

Wave  Function 

1)+    -l)l/x/2 

B 

(-^-1) 

1>-  1  -l>/v/2 

C 

(1,0) 

1  0) 

D 

(L-i) 

1)-     -1)-V2    0)/2 

E 

ilir 

i{   l)+     -l))-\/2    0)/2 

F 

(0,0) 

(-1  +  0     l)  +  (l  +  0     -l)-\/2    0)l/v/6 

G 

i-ho) 

(-1  +  0    i)  +  (i  +  0    -i)/2 

Rx{^)-  i'C  -^  i'A^'4^D  -^  4>G\^B,'^F,-^E  fixed. 
Rz{^)  :  V'A  -*  0/7,  *i,'  -^  V'£>;V'C,0F,V'G  fixed. 


The  i2x(x)  transforms  the  points  {a',ri)  into  new  points  [o" ,r)')  given  by 


1 


1  3 

2  4 
1 


(4.28) 


which  corresponds  to  the  mapping  given  by  eqs.(4.25). 

There  is  some  confusion  in  the  literature  concerning  the  values  of  o' .  In  Van 
Kranendonk's(75)  book  the  negative  value  of  o'  had  been  ruled  out  as  it  was 
unphysical.  But  this  is  not  strictly  correct.  Van  Kranendonk's  remarks  refer 
to  considerations  of  the  pure  states  |  J^  =  0)  and  |  J^  =  ±l)  only.  He  does  not 
discuss  either  the  density  matrix  approach  or  the  formulation  of  the  intrinsic 
quadrupolar  order  parameters  o'  and  r)  needed  to  describe  the  orientational 
degrees  of  freedom  of  the  ortho-i/2  molecules. 

The  results  of  the  relabeling  transformation  are  particularly  easy  to  under- 
stand if  we  consider  the  pair  (5  =  cr',N  =  Y'^)  which  transform  orthonor- 
mally  when  the  axes  are  rotated.    The  parameter  space  {S,  N)  is  shown  in 


fi 


[   -' 


# 
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f37ri 


Fig.13.  The  rotation  ii!a;('^j'')''"  i^^^)  parameter  space  through  the  line  BZ 
with  the  transformed  points  given  by        ■  :        . 


N')  =  {2    f  lU)  ^         ^29) 


S     \         I       2         2      1   I    5 

/3         1 
2  2 

Similarly,  the  relabeling  {x,  y,  z)  ->  (y,  a:,  2)  corresponds  to  a  reflection  through 

the  line  CG.  The  entire  area  of  allowed  values  in  parameter  space  can  therefore 

be  mapped  out  starthig  only  with  the  triangle  ACFE  by  simply  relabelling  the 

principal  axes.  We  only  need  to  consider  the  hatched  region  of  parameter  space 

shown  in  Fig.] 2  in  order  to  describe  all  physically  distinguishable  orientational 

states  for  ortho-if2  molecules  in  the  solid  state.  This  is  not  the  only  choice  that 

can  be  made  for  a  "primitive"  area  of  inequivalent  values  of  {a',ri).   We  may 

also  choose  the  sum  of  ACFK  and  AGFJ  in  Fig.13.  In  that  region  the  states 

have  the  minimum  values  of  the  eccentricity  N.  Both  choices  are  equivalent. 

A  Proposition  for  A  Zero  Field  Experiment 
An  equivalent  expression  for  the  spectrum  for  each  molecule  is 


Ai/i  =  ±y-D(7r2o(i))oz  ':  (4.30) 

As  previously  discussed,  considering  the  transformation  to  the  local  frame 
we  obtained 

Au^  =  ±D[-a'-p2{cos9,)  +  -r)isin'^0,cos2<j)^]  (4.31) 

where  (<?»,  <^i)  are  the  polar  angles  defining  the  orientation  of  the  applied  mag- 


i  netic  field  with  respect  to  the  local  molecular  symmetry  axes. 
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Figure  12.  The  Allowed  Values  of  a  And  r? 
for  Spin-1  Particle 
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Figure  13.  Diagram  of  Allowed  S  And  N 
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We  propose  that  the  assumption  of  axial  symmetry  can  be  tested  exper- 
imentally by  examining  the  zero  field  NMR  absorption  spectrum.  Reif  and 
Purcell  (76)  have  carried  out  zero  field  studies  for  the  long  range  ordered  phase 
where  it  is  known  that  o'  is  constant  and  /?  =  0,  but  it  has  not  previously  been 
considered  for  the  glass  phase. 

In  zero  applied  magnetic  field  the  degeneracy  of  the  nuclear  spin  levels  is 
lifted  only  by  the  intramolecular  spin-spin  interactions. 

HDD{i)  =  hDl-^iSll  -  If)  +  \m{ll^  +  Il^)]  (4.32) 

The  tensorial  operators  for  the  rotational  degrees  of  freedom  have  been  re- 
placed by  the  expectation  values  a'  and  rj.  An  applied  radio-frequency  field 
can  induce  magnetic  dipole  transitions  between  the  nuclear  spin  levels  in  anal- 
ogy with  the  so-called  pure  quadrupole  resonance  absorption (51  Chapter  VIII). 
The  eigenstates  of  H^j)  are  |  /^  =  0)  and  |  ±)  =  ^(|  /^  =  1)±  |  /^  =  —  1). 
For  each  molecule  i,  three  resonance  lines  can  be  expected  corresponding  to 
the  transitions  |  +)  — >|  0),|  — )  — >|  0)  and  |  +)  ^-1  — )  with  frequencies 
I//  =  D{-o'-  +  lr)i),  uf  =  D{-o[  -  \T]i)  and  uj  =  Dr^i,  respectively.  If 
axial  symmetry  is  a  good  approximation,  there  is  only  one  line  at  i/j  =  ~D^i 
and  the  detailed  shape  of  the  NMR  absorption  spectrum  in  zero  field  will  be 
identical  with  that  observed  at  high  fields.  Otherwise  the  high  and  zero  field 
spectra  will  not  be  identical. 

A  Model  for  The  Distribution  Function  of  o 
X.  Li,  H.  Meyer  and  A.  J.  Berlinsky  (23)  have  proposed  a  model  for  the 
distribution  function  of  a'.  They  assume  that  for  a  single  crystal  the  orientation 
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of  the  principal  axes  are  uniformly  distributed.  They  considered  a  Cartesian 
basis  for  the  description  of  the  single  particle  density  matrices 

The  pure  states  p  ■      =  1  correspond  to  ^Q-wave  functions  for  a  =  x,y,  or  z. 

(a) 

This  leads  to  a  natural  parameterization  of  the  p       given  by 

.  "^  - - 


where  the  "energies"  Eai  are  in  general  temperature  dependent.  Li  et  al  (23) 
then  made  the  further  assumption  that  the  effective  site  energies  are  normally 
distributed  about  zero  with  a  width  A(r)  subject  to  the  constraint  J2a  -^"«  ~ 
0.  i.e. 

P{E.,E,,E.)  =  -^,e-m^^l^El)l^^6[E,-^Ey  +  E.) 

These  assumptions  lead  to  definite  predictions  concerning  the  variation  of  the 
NMR  line  shape  parameters  M2  and  M4,  and  in  particular  of  the  variation  of 
-ttI  as  a  function  of  the  degree  of  local  orientational  order  parameter  measured 

Instead  of  Cartesian  symmetry  for  the  effective  site  energies,  we  explored 
the  same  trends  for  a  cylindrical  symmetry  for  the  effective  site  energies  because 
of  the  axial  symmetry  of  the  quadrupolar  interactions.  This  also  leads  to  a 
natural  description  of  the  energy  states  in  terms  of  local  "two  level  system," 
corresponding  to  Jz^  =  ±1  and  Jz^  =  0  separeted  by  an  energy  gap  A^. 
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The  probability  distribution  of  A  is  Gaussian  as  following 


The  order  parameter  a  =(3J^  —  2)  will  be 


^(A)  =  \/^e    2^r,'2  (4.33) 


(a(A))  =  '^"t      '  ,  (4.34) 

2e    KT  +1 


We  can  prove  that 

•  M4        15    (a^) 


Mi  7   ((^2))2 


(4.35) 


where 


let  §  =  X 


cr2    =   /       ^^ LJ  ^-A   ^^  use) 

,  o,         f°°  (2e~Kf^  -2)^    [2    -x'^  ,        , 

^0      (2e~K^^  +  1)4  »  ^ 
By  using  numerical  integral  method  the  curve  in  Fig.  14  shows  the  ratio 

'     '2  versus  \/(cr2) (order). 

Considering  the  intermolecular  and  the  intramolecular  dipolar  broadening, 

the  fourth  moment  M4  and  the  second  mement  M2  should  be  (51) 

-^  -  -'-W 


where 


yl  ytntra   .    ymter 

^1  =A^l"'l^  3 \^hh-I     ■I) 

k,l  ^kl 


61 

;/."'-  =  7..,*  E—^^^.si 

i.j  ij 

Following  a  quite  complicated  but  straightforward  calculation  we  can  prove 

that 

M2  =  Mi""^'"  +  A^"*'"  (4.41) 

M4  =  Mf  *'■«  +  Ml""^''  +  4M^"*'''^m|"*^'"  (4.42) 

For  a  Gaussian  shape  one  has  Ml"^^*"  =  'i{M'rP^^^)'^  and 

M4  =  MX'''''  +  3(M^'^*^'')^  +  4i^'^*'''^M^"*^'"  (4.43) 

Taking  the  same  value  of  M^'***'"  =  20(A;/i2;)^  as  in  ref.23  we  obtain  the  curve 

\jintra  . — —  w 

in  Fig. 15.  The  upper  one  is  fj^tntra\2  versus  y  {o^)  and  the  lower  one  is  tjt^ 
versus  y^o"^). 

The  experimental  lineshapes  (23,24,27)  are  below  the  calculated  lineshape, 
but  the  shape  and  slope  are  almost  the  same.  The  value  of  maximum  -^  of 
experimental  lineshape  is  3.15  (18)  and  the  value  of  maximum  -772  of  calculated 
lineshape  is  5.60. 
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CHAPTER  5 

NMR  PULSE  STUDIES  OF  SOLID  HYDROGEN 
The  experimental  technique  of  nuclear  magnetic  resonance  with  spin  echoes 
has  been  used  widely  in  recent  years  as  a  tool  to  investigate  the  dynamical 
properties  of  crystals.  In  particular,  much  attention  has  been  paid  to  the 
study  of  relaxation  rates  of  dynamical  processes  in  molecular  solids  by  means 
of  the  spin  echoes  occuring  after  application  of  resonant  rf  90°  —  r  —  V'U)  pulse 
sequences  and  rf  90°  -  r  -  '4)i^\  -  tw  -  ^j't.n  pulse  sequences. 

The  purpose  of  this  chapter  is  to  develop  the  theory  of  spin  echoes  of  solid 
H2  in  orientationally  ordered  phase  and  discuss  the  low-frequency  dynamics  of 
orientational  glasses. 

Solid  Echoes  ' 

There  are  several  publications  which  have  presented  calculations  for  the 
amplitude  of  solid  echo  responses  to  90°  —  t  ~  ip^  pulse  sequences  (39,  78,  79). 

We  would  like  to  examine  the  time  dependence  of  the  nuclear  operators 
and  the  conditions  for  the  focussing  of  solid  echoes. 
The  Time  Dependence  of  The  Nuclear  Operators 

A  general  two-pulse  sequence  denoted  by  (f  - 1"  -  ipcf,)  is  sketched  in  Fig. 16. 

In  NMR  experiments  one  observes  only  the  ortho-molecules  (7  =  1).  In  a 
first  approximation  we  will  treat  the  system  as  a  set  of  independent  molecules. 
The  signal  we  measured  is  proportional  to  Tr[p{t]I+]  =  Tr{p{t){Ix  +  Hy)], 
where  p  is  the  density  operator  at  time  t  and  /  the  nuclear  spin  operator. 
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Figure  16.  Solid  Echo  (two-pulse  sequence) 
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Figure  17.  Stimulated  Echo  (three-pulse  sequence) 
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In  equilibrium  the  expression  of  density  matrix  is  given  by 


(5.1) 


Let  Z  be  the  direction  of  applied  magnetic  field. 

where  lOq  is  the  Larmor  frequency. 

Since  ^  <  1,  e^^"  «  1  +  ^hi 

P^e,  «     ^'  ''  %^"^     «  -(1  +  ^/,,)  (5.2) 

""'       Tr{l  +  ^I,,)       S^        KT  '''  '      ' 

We  know  that  the  radio  frequency  field  is  responsible  for  nonequilibrium 

behavior  of  the  system. 

In  the  rotating  frame  after  the  first  90°  pulse,  the  density  matrix  p{0+) 

would  be 

p(0+)  =  e-UH^'^*p^{0_)eUH^dt  (5  3j 

If  the  radio  frequency  field  is  in  the  y  direction  of  the  rotation  frame,  we 
obtained 

M0+)  =  eJ5'.|l(l  +  |^y|.-JI'. 

p{0+)  =  ^(1  +  ^/.)  (5.4) 

This  represents  that  the  first  pulse  Rq  puts  the  magnetization  along  the  x  axis 
(for  simplicity  we  have  dropped  the  site  index  i). 

We  shall  consider  the  time  evolution  of  the  nuclear  spin  operators  under 
the  effect  of  the  intramolecular  Hamiltonian  in  the  rotating  frame. 

Hdd   =  E  -f'^^Iz^+  E  \\llDhP2{0,)ajf,  (5.5) 
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where 

4  =  ^^(34-2) 

and 


3  KT 

It  is  assumed  that  the  second  pulse  is  a  V'—  pulse  and  the  pha^e  of  that  is 
(j)  (the  angle  between  radio  frequency  field  and  y  axis  in  rotating  frame).  After 
the  second  pulse  the  density  matrix  is  given  by 

p(t+)  =  e'^i'{^=^^'^n't>+IyCos4>)    t_\-irl){I:,sin<j>+IyCos<l>)  f^y\ 

The  nuclear  spin  density  matrix  at  time  t  after  second  pulse  is  therefore 
{t  +  T)  =  u{t)R{.^)u{T)l{l  +  ^/.)u^(r)it:t(<^)«t(f)  (5.8) 


Pk-     .     •/  -v-,--.T.-v    .3V-         J^J, 


where 


and  the  signal 


R{(f>)  =  e^W-'>^n4>+^v<'os4>)  (5  10) 

S{t  +  T)  (xTr{p{t  +  T)I+} 


Tr{p[i  +  T)I^} 

1  hbJi 


'^Tr{u{t)R{<l>)u{T)hu\T)R\(l>)u\t)I+} 


3KT 

=  3^     X^    (mi  |u(0  |m2){m2  |i?(<^)  |m3)(m3  |u(r)  |m4) 
mi— ms 

(m4  I  Ix  I  m5)(m5  |  u\t)  \  mQ){mQ  \  R^{<p)  \  m-j) 
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{m-j  I  u'(i)  I  mg,){m^  \  I-^  \  mi)  (5.12) 

Formation  of  Solid  Echoes 

After  a  lengthy  calculation  the  results  of  eq.(5.12)  are  still  quite  compli- 
cated (for  general  (j)  and  ^) .  For  some  special  values  of  (f)  and  tp  the  results  are 
the  same  as  that  in  ref.39.  We  would  like  to  give  the  results  in  some  special 
conditions. 

{l)il)  =  0  (apply  one  pulse) 

Tr{p{t)I+}  =  l^e-^'"^'cos[^-DP2{6)at]  (5.13) 

This  is  called  the  Free  Induction  Decay. 

(2)  </>  =  CV^  =  f  This  is  a  90°  -  r  -  90°  pulse  sequence 

Tr{p{t  +  T)I+} 

=  l^e-^^-^cos[lDP,{e)a{t  -  r)](l  -  e^^^^rj  (^  ^^j 

Obviously,  Tr{p{t  +  t)I+}  has  a  maximum  value  at  t  =  r.  This  is  what 
is  the  meant  by  a  solid  echo.  Another  question  is  that  according  to  eq.(5.14) 
the  in-phase  echo  (Aw  =  0)  does  not  exist.  This  is  in  disagreement  with  the 
experimental  results.  The  reason  lies  in  neglecting  intermolecular  dipole-dipole 
interaction  and  the  difficulty  of  getting  exact  Aw  =  0  in  experiments. 

(3)<^  =  f ,  V  =  f  This  is  a  90°  -  7-  -  90°  pulse  sequence 

Tr{p{t  +  r)I+} 
=  l^e-^''^'cos[\DP2i0)ait  -  r)](l  +  e^'^^^rj  ^^^^^ 

From  eq.(5.15),  same  as  the  second  case,  Tr{p{t  +  T)I^}  obtaines  maximum 
value  at  t  =  T,  forming  the  solid  echo. 
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Stimulated  Echoes 

It  is  already  shown  that  in  certain  cases  a  sequence  of  three  90°  pulses  may 
be  advantageous  (80,81). 

We  first  describe  the  formation  of  nuclear  spin  stimulated  echoes.  The 
stimulated  echoes  can  be  used  to  compare  a  "fingerprint"  of  the  local  molecular 
orientations  at  a  given  time  with  those  at  some  later  time  (less  than  Ti)  and 
thereby  used  to  detect  ultra-slow  molecular  re-orientations. 

Applications  to  the  study  of  the  molecular  dynamics  on  cooling  into  the 
quadrupolar  glass  phase  of  solid  hydrogen  will  be  discussed. 
Formation  of  Stimulated  Echoes 

A  three-pulse  sequence  is  sketched  in  Fig. 17.  As  the  analysis  of  solid  echoes, 
the  nuclear  spin  density  matrix  at  different  time  are  given  by 

;,(0^)  =  p.,  «  i(l  +  ^h) 

/9(r_)  =  e    ''•'0     'id  y  I     ;o(0-j-)e^  •'o     <*<' ^  ' 
p{t+)  =  ^iWxsin4>-\-IyCos<i))    i_\-i\l>{I^sin4,-\-IyCos(j>) 

p[{tw  +  r)-]  =  e    h"dd^^p[Tjf.)eh"dd^«> 
p[[tw  +  r)+]  =  e''^'(^-"*"^'+'^'''^')^[(<^  +  r)_]e-»V''(4sm^'cos<^') 


p{tw  +  r  +  t)  =  e    h"'^d'-p[{tu)  +  T)  +  \eh 


dd 


The  signal 

S{tu>  +  T  +  t)<xTr{p{tu,  +  T  +  t)I+}  (5.16) 

The  first  preparatory  pulse  creates  a  transverse  magnetization  in  the  rotat- 
ing frame.  Under  the  influence  of  the  Hamiltonian  given  by  equation(5.15),  the 
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evolution  during  the  short  time  r(<C  T2)  leads  to  the  formation  of  nuclear  spin 
states  described  by  p{t-)  which  contains  both  transverse  magnetization  and 
transverse  alignment.  The  transverse  components  are  transferred  by  a  second 
pulse  i2(V')  into  longitudinal  components  corresponding  to  spin  polarization 
and  spin  alignment.  They  will  be  stored  and  evolved  during  long  waiting  time 
tw  (chosen  short  compared  with  the  longitudinal  relaxation  times).  Therefore 
we  can  obtain  a  "fingerprint"  of  the  local  alignments  o.  After  the  waiting 
period  iu;  ^  T2  the  stored  components  can  now  be  'read'  by  means  of  a  third 
pulse  of  rotation  angle  0'.  The  spin  polarization  and  spin  alignment  'stored' 
during  tw  are  transferred  by  the  third  pulse  into  coherent  transverse  states, 
which  then  evolve  in  a  reverse  manner  to  that  occurring  during  the  first  evo- 
lutionary period  r  and  the  transverse  signal  focuses  to  stimulated  echo  after  a 
delay  time  r  following  the  third  pulse. 

Tr{p{tuj  +  T  +  t)I+}  ! 

=  -^     2^     (mi  I  e    htiad     (  m2)(m2  I  R{(f>')  \  mz){mz  \  e    h^dd  «-  |  m^) 
mi— mi2 

(m4  I  R[<t>)  I  m^){m^  \  e~h^dd'^  \  mQ){mQ  \  Ix  \  m-j) 

{mj  I  eh^cui  ^  I  ms){ms  \  R{4>)  \  mg){mQ  \  eh"dd  ^^  |  mio) 

(mio  I  i?+(<^')  I  mii)(mii  |  e^^'"  *  |  mi2)(mi2  |  /+  |  mi)  (5.17) 

Since  the  calculation  is  extremely  tedious  and  the  results  for  general  (^,0' 
and  general  V'?V''  are  very  complicated,  we  only  investigate  the  results  of  cal- 
culations by  taking  two  special  cases. 

(1) 
0  =  0,  <^'  =  f 

0  =  f ,  V'  =  f 
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This  is  a  90°  —  r  —  90°  —  t^j  —  90°  pulse  sequence 

Tr{p{t^  +  T  +  t)I+} 

I 

=  l^e-^^''Ucos{Au;T)cos^Au;tu^)cos{lDP2{e)a{t  -  r)] 
Z  lil  z 

+  --^e~"'^'^Ucos{AujT)sin^{Aujtu))cos\-DP2{0)a{t  +  t)] 
2  KT  2 

sin{AujT)sin{Au)tw)cos[-DP2{0)a{t  -  t^  +  r)]  (5.18) 


—  ■ 

2  hbjQ  _. 
3KT^ 

iAwt 

(2) 

^ 

=  h^' 

_  T 
-2 

0 

=  h^' 

5r 
-  2 

This  is  a  90°  —  r  —  90^  —  ^ly  —  90°  pulse  sequence 
Tr{p{tr^  +  T+t)I+} 

=  -^e-'^^Hsin^[Autu,)sin[AuT)cos\-DP2[e)o[t  -  t)\ 

+  ^|^e~*^'^*co5(AwOcos(Au;i«;)co5[jDF2W^(f  -tw^  r)\  (5.19) 

Examining  expressions  of  eqs.(5.18)  and  (5.19),  besides  the  echo  at  i  =  r 
after  the  third  pulse,  there  is  a  additional  echo  at  t  =  tw  —  t.  The  results 
are  sketched  in  Fig. 18.  The  appearance  of  multipole  echoes  has  been  seen  in 
experiments  (Fig.  19)  for  ortho-para  hydrogen  mixtures  at  low  temperatures. 

Engelsberg  et  al.  (81)  have  presented  results  for  nuclear  spin  stimulated 
echoes  in  glasses.  The  curve  of  ^^B  echoes  in  borosilicate  glass  at  4.2  K  (Fig. 20) 
shows  that  in  addition  to  the  solid  echo  at  f  =  2r (after  first  pulse)  and  a 
stimulated  echo  at  t  =  T  +  2t,  other  echoes  at  t  =  2T  (image  echo)  and 
t  =  2T  +  2r  (primary  echo)  were  clearly  observed.  For  longer  waiting  times, 
the  solid  echo  (which  they  called  spontaneous  echo)  decays  rapidly  and  only  the 
stimulated  echo  remains  detectable. 
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Fig. 21  shows  the  temperature  dependence  of  the  (f  —  t"  —  f  —  ^ly  —  f ) 
stimulated  in  the  quadrupolar  glass  phase  of  solid  hydrogen  (r  =  25fj,s,  tw  = 
2ms).  The  theoretical  results  (eqs  5.18  and  5.19)  and  the  experimental  results 
clearly  show  that  the  amplitude  of  stimulated  echo  is  proportional  to  jj.  For 
experimental  curve  only  a  very  slight  modificaion  (indicated  by  the  arrow)  was 
observed.  We  will  discuss  this  phenomenon  later. 

The  experimental  curve  in  Fig. 22  gives  the  relation  of  stimulated  echo 
versus  waiting  time  tu)  of  solid  hydrogen  (ortho-concentration  x  =  0.54,  T  = 
220mK).  The  time  scale  is  logarithmic.  The  logarithmic  decay  behavior  can 
be  understood  in  terms  of  the  motional  damping  of  stimulated  echoes. 

Low-Frequency  Dynamics  of  Orientational  Glasses 
The  orientational  glasses  (20,  solid  ortho-para  H2  mixtures  (6),  N2/A  mix- 
tures (83)  and  the  KBri_j.K{CN)x  mixed  crystals  (82,  84,  85))  form  a  sub- 
group of  the  general  family  of  spin-glasses  which  continue  to  generate  intense 
interest  because  of  the  apparently  universal  low  temperature  properties  ob- 
served for  a  very  diverse  range  of  examples  (dilute  magnetic  alloys,  mixed  crys- 
tals, dilute  mixtures  of  rotors,  partially  doped  semiconductors  (86),  Josephson 
junction  arrays  (87)  and  others).  The  most  apparent  striking  universal  fea- 
tures (20)  are  an  apparent  freezing  of  the  local  degrees  of  freedom  on  long 
time  scales  without  any  average  periodic  long  range  order,  characteristic  slow 
relaxations  and  history-dependence  following  external  field  (magnetic,  electric, 
elastic-strain....)  perturbations,  and  a  very  large  number  of  stable  low  energy 
states. 
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Figure  18.  Sketch  of  The  Results  of  Calculations 
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Figure  19.  Experimental  Curve  (z  =  23%,T  =  38mK) 
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Nuclear  spin  stimulated  echoes  in  glasses 
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Figure  20.  ^^B  Spin  Echoes  in  Borosilicate  Glass  at  4.2K 
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Figure  22.  The  Observed  Decay  of  Stimulated  Echo 

(x  =  0.54,  T  =  220mK) 
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The  echo  calculation  mentioned  above  was  based  on  the  static  case.  If  the 
ortho  molecules  are  in  slow  motion,  the  stimulated  echoes  will  be  damped.  At 
low  temperatures  the  random  occupation  of  lattice  sites  for  solid  H2  mixtures 
(for  X  <  55%)  leads  automatically  to  the  existence  of  local  electric  field  gra- 
dients, the  field  conjugate  to  the  local  order  parameter,  which  plays  the  same 
role  as  the  magnetic  field  for  the  dipolar  spin  glasses.  This  random  local  field 
therefore  makes  the  problem  of  local  orientational  ordering  in  random  mix- 
tures equivalent  to  the  local  dipolar  ordering  in  spin  glasses  in  the  presence  of 
random  magnetic  fields.  i 

In  analogy  with  the  analysis  for  spin  glasses  we  assume  that  the  existence  of 
local  electric  field  gradients  leads  to  clusters  (or  droplets)  of  spins  (88).  Based 
on  Fisher  and  Huse's  recent  picture  (89),  we  provide  a  explanation  of  the  low- 
frequency  relaxation  and  the  low-temperature  specific  heat  of  solid  ortho-  para 
hydrogen  mixtures. 

In  the  scaling  model  of  Fisher  and  Huse  the  low-energy  excitations  which 
dominate  the  long-distance  and  long-time  correlations  are  given  by  clusters  of 
coherently  reoriented  spins.  Their  basic  assumptions  are: 

(1)  Density  of  states  at  zero  energy  for  droplets  (d  dimension)  length  scale 
L  as  L-\  where  0  <  ^  <  ^. 

(2)  Free  energy  barriers  Eq  for  cluster  formation  scale  as  Eg  ^  L^  with 
e  <rp  <d-l. 

With  these  assumptions,  Fisher  and  Huse  show  that  the  autocorrelation 
function 

decays  as  {\ogt)    "Z"  for  f  ^  cx). 
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For  our  system,  assuming  axial  symmetry,  the  quasi-static  local  orienta- 
tional  order  parameters  are  the  alignments  o^  =(3  j|  —  2)^-  and  the  correspond- 
ing autocorrelation  functions  Cj(^)  =  Oj^[Qi)o^[i)  can  be  studied  directly  by 
NMR. 

For  a  90°y  —  r  —  90°j/  —  t^  —  90°  pulse  sequence  we  assume  at  <  =  r,  order 
parameter  a  —  a{T)  and  at  ^  =  r  +  tw,o  =  o{t  +  iw^  for  each  ortho  molecule 

3  i\  J  2  2 

3    I\  1  id  il 

3  J\.i 

cos^-DP2{Q)o{t  +  fu,)i  -  \dP2{Q)o{t  +  «^)^u,  +  ^£>P2 W^(0^1        (5-21) 
2t  £i  ^ 

Considering  <  =  r,  the  stimulated  echo  amplitude 

A  oi{{cos\D,To,[r)\cos\D,Ta,{T  +  t^)\))  (5.22) 

where  the  double  brackets  refer  to  an  average  of  configuration  and  D^  = 
\DP2{cosB^.  The  important  point  is  that  if  the  local  order  parameters  a^ 
remain  fixed  during  txu-,  there  is  no  damping  of  the  stimulated  echo,  while 
0{  changes  due  to  local  re-orientations,  then  the  contribution  to  the  echo  is 
severely  attenuated.  The  product  Dr  can  in  practice  be  made  very  large  and 
this  method  can  therefore  be  used  to  study  ultra-slow  motions  in  solids.  We 
believe  that  in  Fig. 21  the  departure  portion  from  j.  (indicated  by  an  arrow)  is 
due  to  slow  motion. 

A  barrier  will  have  a  characteristic  life  time  given  by  an  Arrhenius  Law 


1         1     ^ 


—   =   t^B"^ 
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or  tunneling  rate  T 


T{EB)  =  Toe    ^B-r  (5.23) 

where  Tq  is  the  characteristic  attempt  frequency  for  clusters  of  this  size.  In 
the  long  time  limit  Tq  is  resonably  well-defined  because  it  is  associated  with  a 
characteristic  cluster  size.  In  a  time  t  the  only  barriers  crossed  will  be  those  sat- 
isfying 0  <  Eb  <  Emax{t)  where  Emax{t)  =  K^T log  f.  Any  barriers  crossed 
lead  to  significant  changes  in  the  local  order  parameters  and  the  amplitude  of 
the  stimulated  echo  is  then  simply 

/•oo 

A{t)  =   /  P{EB)dEB  (5.24) 

At  low  temperatures,  assuming  a  constant  density  of  barrier  heights  P{Eb), 
we  find  I 

A{t)  =  l-KBTPo\og{^)  (5.25) 

The  prefactor  Pq  can  be  determined  from  the  low  temperature  behavior  of  the 
heat  capacity. 

For  the  ortho-/f2  molecules  with  angular  momentum  J  =  1,  we  can  asso- 
ciate a  simple  two  level  system  with  the  energy  states  for  a  given  molecules; 
the  states  J^.  =  ±1  being  separated  from  the  state  J^.  =  0  by  a  gap  3Aj  (The 
states  Jf-  =  ±1  are  degenerate  if  there  are  no  interactions  which  break  time 
reversal  symmetry).  At  low  energies  we  can,  following  the  above  arguments, 
identify  the  low  energy  excitations  (which  determine  Cv  at  low  T]  with  a  broad 
quasi-constant  distribution  P(A)  for  0  <  A  <  Aq  for  the  spins  in  a  cluster. 
Identifying  P(A  =  0)  with  Pq,  the  density  of  low  energy  barriers,  we  have 

Cv         „     f^°     18  a2 


^0 


NxR  ju      ^^"^    4e    KbT  ^gK^T  ^4 


where  x  is  the  ortho-i/^2  concentration. 
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let  u  = 

-  Q     A,, 

'"^KbT 

Cv     _       2KbT    f^K^r             u2 
NxR  "    "^      3      Jo           4e-«  +  e«  +  4 

set  t  = 

■W 

NxR         ""    ^3"  Jo    4e-«  +  e«  +  4 

let 

,       2     fT          v}du 

^       3   Jo    4e-"  +  e«  +  4 

du  (5.27) 


(5.28) 


(5.29) 

The  resulting  C'^  (Figure  23)  has  a  linear  temperature  dependence  at  low  T 
and  a  peak  at  T„)^  =  0.70-^  in  close  resemblance  to  the  temperature  behavior 
observed  by  Haase  et  al.  (90).  From  the  peak  position  in  the  experimental 
data,  ^  =  1.27,  PqK  =  0.86  and  for  the  stimulated  echo  decay  this  value 
gives 

Aa/c.(0  =  l-0.43logio(f)  (5.30) 

and  the  observed  decay 

Aobs.it)  =  l-0.55\ogio{r)  (^-31) 

at  r  =  0.22K  for  x  —  54%.  The  agreement  is  remarkably  good. 

The  experimental  curves  indicate  Iq  ~  10~^s. 

It  should  be  noted  that  the  argument  relating  the  logarithmic  decay  to 
maximum  barrier  height  crossed  in  time  t  can  also  apply  (over  a  short  time 
scale)  to  the  case  of  orientational  ordering  in  pure  N2  studied  in  reference 
(11,91)  because  one  also  observes  a  relatively  large  distribution  of  order  pa- 
rameters centered  on  a'  =  0.86  and  with  width  0.12  in  this  case.  The  essential 
point  is  that  the  time  scale  of  the  slow  relaxations  in  the  glass  phase  is  simply 
related  to  the  low  temperature  behavior  of  the  heat  capacity. 
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Another  important  point  is  that  the  characteristic  times  Iq  are  much  shorter 
than  the  spectral-diffusion  time  scale  (~  sec)  seen  by  the  recovery  of  holes 
burnt  in  the  NMR  lineshape.  We  therefore  find  it  difficult  to  attribute  the 
logarithmic  decay  seen  in  H2  to  spectral  diffusion  across  the  NMR  linshape. 
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Figure  23.  Calculated  Curve  oi  C'^  -  t 


CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 
The  ortho-para  mixtures  of  solid  H2  are  studied  theoretically  for  fee  lattices 
(X  >  0.55)  of  finite  site  by  Akira  Mishima  and  Hiroshi  Miyagi  (92).  The  sys- 
tematic theoretical  studies  of  nuclear  magnetism  for  quadrupolar  glass  regime 
are  carried  out  in  this  dissertation. 

We  have  developed  a  theory  of  the  nuclear  spin-lattice  relaxation  of  ori- 
entationally  ordered  ortho  hydrogen  molecules  for  the  case  of  local  ordering 
in  the  quadrupolar  glass  phase  of  solid  hydrogen.  We  have  investigated  the 
temperature  dependence  of  T^.  It  shows  that  Gaussian  Free  Induction  Decay 
is  a  quite  good  approximation.  The  calculations  indicate  a  strong  spectral  in- 
homogeneity  of  the  relaxation  rate  rj^^(Ai/)  throughout  the  NMR  absorption 
line.  The  detailed  dependence  is  much  stronger  than  the  simple  dependence 
rr^(Ai/)  oc  (2  +  ct)  given  by  earlier  estimates  by  A.  B.  Harris  et  al.  (1).  There 
is  also  a  strong  variation  with  the  orientation  of  the  local  symmetry  axis  with 
respect  to  the  external  field.  The  variation  is  in  agreement  with  that  found 
by  Hardy  and  Berlinsky  (93)  for  the  long-range  ordered  Pa^  phase.  If  the 
cross-section  relaxation  between  different  isochromats  is  taken  into  account, 
the  results  of  calculations  of  spectral  inhomogeneity  is  in  good  agreement  with 
the  experimental  results. 

In  addition  to  the  spectral  inhomogeneity,  the  relaxation  is  also  found 
to  be  strongly  nonexponential.  This  can  be  understood  easily.  As  a  result 
of  the  glassy  nature  of  the  system,  there  is  a  broad  distribution  of  both  local 
orientational  order  parameters (cr)  and  the  orientation  (a)  of  the  local  symmetry 
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axes  for  the  molecular  alignments.  The  existence  of  these  distributions  at 
low  temperatures  means  that  a  given  frequency  Aw  in  the  spectrum  comes  from 
all  the  different  allowed  combinations  of  a  and  a  that  satisfy  Au  =  2DP2{oc)o. 
The  dependence  of  the  relaxation  rates  on  a  and  P2  then  leads  to  a  distribution 
of  local  relaxation  rates  for  a  fixed  Au.  It  is  this  distribution  in  rates  that  leads 
to  the  observed  non-exponential  behavior  of  M{t).  Calculations  based  on  the 
expected  probability  distributions  for  o  and  P2  yield  results  in  good  agreement 
with  the  results  reported  in  the  publications  (14). 

The  most  important  conclusion  of  the  study  of  Ti  is  that  the  spectral  inho- 
mogeneity  and  nonexponential  recovery  both  results  in  an  order  of  magnitude 
variation  of  the  nuclear  spin  relaxation  and  must  therefore  be  correctly  under- 
stood and  accounted  for  before  attempting  to  analyze  the  experimental  results 
in  terms  of  the  fundamental  molecular  motions. 

We  have  determined  that  in  general  the  orientational  degrees  of  ortho- 
H2  molecules  in  the  solid  need  to  be  described  in  terms  of  density  matrices. 
The  ortho  molecules  have  momentum  J  =  1  and  the  single  particle  density 
matrices  are  completely  determined  by  five  independent  parameters  (if  the 
angular  momentum  is  quenched).  These  parameters  are 

(1)  the  three  principal  axes  (x,y,z)  for  the  second  order  tensor. 

(2)  the  alignment  o'  =(1  -  |  J|),  and 

(3)  the  eccentricity  rj  —{Jx  ~  Jy)- 

The  positivity  conditions  for  the  density  matrix  show  that  the  only  allowed 
values  of  (a',??)  are  those  enclosed  in  a  triangle  in  [o' ,r))  space  whose  vertices 
are  the  pure  states  |  Jz  =  0)  and  |  J^  =  ±1).  Not  all  of  these  allowed  values 
are  physically  inequivalent  because  one  may  relabel  the  principal  axes  and  we 
have  shown  that  one  can  determine  a  simple  primitive  set  of  order  parameters 
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which  are  inequivalent  by  the  choice  2o'  >  r]  >  0.    Orientational  states  with 
negative  o'  are  not  excluded  on  theoretical  grounds. 

Studies  of  molecular  dynamics  have  shown  the  existence  of  para-librons, 
i.e.  collective  excitations  that  exist  in  large  clusters  with  well-developed  short 
range  order. 

Nuclear  spin  stimulated  echoes  have  proved  to  be  very  effective  for  the 
study  of  ultra-slow  molecular  motions  both  in  the  molecular  orientational 
glasses  and  in  ordinary  glasses  (81).  The  existence  of  multiple  echoes  which  is 
given  by  theoretical  calculations  is  also  in  good  agreement  with  experimental 
results  of  solid  hydrogen  and  ordinary  glass  at  low  temperature. 

We  have  offered  a  unified  explanation  of  the  slow  relaxational  behavior  and 
low  temperature  heat  capacity  of  the  quadrupolar  glass  phase  of  solid  hydrogen 
in  terms  of  the  density  of  low-energy  excitations  in  the  system. 

Since  the  stimulated  echoes  are  damped  by  any  change  in  the  orientational 
states  during  t^,  they  can  be  used  to  detect  very  slow  molecular  motions.  The 
phenomeno  of  logarithmic  decay  of  the  stimulated  echo,  which  is  similar  to 
decay  of  magnetization  in  metal  spin  glass,  can  be  explained  by  using  Fisher 
and  Huse's  recent  picture  of  the  short  range  spin-glasses  and  the  domination 
of  the  long-term  relaxation  by  low  energy  large-scale-cluster  excitations. 

In  analogy  with  spin  glasses  (e.g.  alloys  like  Cu  —  Mn  where  a  random 
configuration  of  spins  condenses  at  low  temperature)  the  linear  specific  heat 
with  temperature  at  low  T  can  be  understood  in  terms  of  the  suggestion  (94) 
that  the  dominant  contribution  to  the  specific  heat  will  be  from  clusters  of 
particles  for  which  the  energy  barrier  is  suflficiently  great  so  that  resonant 
tunneling  between  the  two  local  minima  does  not  occur,  but  sufficiently  small  so 
that  tunneling  between  the  two  levels  can  take  place  and  thermal  equilibration 
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can  occur  during  the  time  span  of  the  specific  heat  experiment.  We  obtained 
a  very  good  linear  dependence  curve  at  low  T  by  using  a  numerical  integral 
method. 

Having  combined  with  the  results  of  logarithmic  decay  of  stimulated  echoes 
and  linear  low-temperature  specific  heat,  a  very  good  numerical  expression  of 
the  amplitude  of  the  stimulated  echo  as  a  function  of  waiting  time  t^  for 
X  =  54%,  T  =  0.22K  is  obtained. 

Most  people  agree  that  molecular  orientations  of  random  ortho-para  hydro- 
gen mixtures  become  frozen  at  low  temperature  and  there  is  no  evidence  for 
any  well-defined  phase  transition  on  cooling  from  the  completely  disordered 
high  temperature  phase.  However,  the  gradual  transition  to  the  glass  state 
in  these  systems  involves  strong  co-operativity  as  evidenced  by  Monte  Carlo 
calculations.  Studies  of  the  nuclear  relaxation  and  molecular  dynamics  have 
shown  quantitatively  how  important  the  co-operativity  is.  Further  theoretical 
study  of  the  slow  motions  needs  to  be  carried  out  to  improve  the  understanding 
of  the  nature  of  the  quadrupole  ordering  in  random  mixtures  at  low  tempera- 
tures and  their  relation  to  the  family  of  glasses.  Further  experiments  at  very 
low  temperatures  will  also  provide  a  deeper  understanding  of  the  properties  of 
the  quadrupolar  glass  state. 


APPENDIX 

FLUCTUATION-DISSIPATION  THEOREM 
The  fluctuation-dissipation  theorem  is  extremely  important  because  it  re- 
lates the  response  matrix  to  the  correlation  matrix  for  equilibrium  fluctuations. 
In  Chapter  2  we  mensioned  the  definitions  of  the  response  function(eq.2.32) 
and  the  relaxation  function(eq.2.35).  The  response  function  my  be  written  in 
a  number  of  equivalent  ways: 

fkl{t)  =  ^Tr{[pe„Mi{-t)]Mk}  (A.l) 

tn. 


^{[Mi{-t),Mk]) 
in 


(A.2) 


fkl{t)  =  ^{[Mi,M,,{t)])  (A.3) 

Using  an  identity  due  to  Kubo  the  response  function  may  also  be  written 
as 


fklit)  =  I    ds{Mi[-ihs)Mk{t))  (A.4) 

Jo 

fklit)  =  -  f    ds{Mi{-ihs)Mk{t))  (A.5) 

Here  ihs  plays  the  role  of  time.  Such  formulas  appear  in  quantum  statistical 
mechanics  because  of  the  formal  similarity  between  the  way  time  and  temper- 
ature must  be  treated.  In  the  calculation  of  a  canonical  partition  function  one 
must  take  into  account  the  factor  e~^" .  In  calculating  the  time  evolution  of 
operators  one  must  consider  operators  of  the  form  e^^^  .By  wntmg  t  =  ins 
one  can  see  that  s  will  play  a  role  identical  to  jS  in  all  formal  developments. 


■:.>{ 


89 


90 

From  eq.(2.35)  and  eq.(A.5)  it  follows  that 

Fki{t)=        ds{Mi{-ihs)Mk{t))  -    \im    /     ds{Mi{-ihs)Mk{T))       (A.6) 
Jo  T-^ooJo 

As  T  apparoches  infinity  in  the  second  term  of  Eq.(A.6)  we  can  assume 
that  the  correlation  is  lost  between  various  components  of  the  magnetization 
so  that  the  correlation  function  factors  into  (M/)(M^).  Hence  the  relaxation 
function  may  be  written  in  its  more  usual  form 

Fkiit)  =  /     ds{AMi{-ihs)AMk{t))  (A.7) 

Jo 

In  this  form  there  exists  a  mathematical  identity  between  the  relaxation  func- 


tion and  the  correlation  function  gki{t),  given  by 

•+00 


/+00 
B{t  -  r)gki{T)dT  (A.8) 

-oo 


where 


j7fci(0  =  ^([AMfcW,AMi(0)]+)  (A.9) 


the  +  subscript  indicates  the  anticommutator;  or  in  terms  of  Fourier  transforms 

hiH  =  %^  (A.io) 

B{u) 


where 


and 


^(^)  =  ^^^^[^^^^(w^i  ^^-''^ 


^  =  ^,,th  (^)  (A.12) 

The  derivation  of  these  identities  is  called  the  quantum-mechanical  fluctuation 
dissipation  theorem.  This  theorem  relates  the  response  of  step  function  distur- 
bance to  the  correlation  function  of  fluctuations  in  the  equilibrium  ensemble. 
In  the  classical  limit  {%  -^  0)  B{t)  ->  l3S{t),  and  one  obtains  the  classical 
fluctuation-dissipation  theorem 

Fklit)  =  I39kl{t)  (A.13) 
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